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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as 
follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

inches 

0.0254 

meters 

pounds  (force) 

4.448222 

newtons 

pounds  (mass)  per  cubic  inch 

27,679.9 

kilograms  per  cubic  centimeter 

pounds  (mass)  per  square  foot 

4.882428 

kilograms  per  square  meter 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

tons  (force)  per  cubic  foot 

32036.9 

kilograms  per  cubic  meter 

tons  (force)  per  square  foot 

95.76052 

kilopascals 
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1  Introduction 


Purpose 


The  purpose  of  this  manual  is  to  provide  a  detailed  discussion  of  techniques 
used  for  the  design/analysis  of  pile  foundations.  Several  of  the  procedures  have 
been  implemented  in  the  CASE  Committee  computer  programs  CAXPILE 
(Dawkins  1984,  Mosher  et  al.  1997),  CPGA  (Hartman,  Jaeger,  Jobst,  and  Martin 
1989)  and  COM624  (Reese  1980).  Theoretical  development  of  these  engineer¬ 
ing  procedures  and  discussions  of  the  limitations  of  each  method  are  presented. 


Pile  Behavior 

The  purpose  of  a  pile  foundation  is  to  transmit  the  loads  of  a  superstructure  to 
the  underlying  soil  while  preventing  excessive  structural  deformations.  The 
capacity  of  the  pile  foundation  is  dependent  on  the  material  and  geometry  of 
each  individual  pile,  the  pile  spacing  (pile  group  effect),  the  strength  and  type  of 
the  surrounding  soil,  the  method  of  pile  installation,  and  the  direction  of  applied 
loading  (axial  tension  or  compression,  lateral  shear  and  moment,  or  combina¬ 
tions).  Except  in  unusual  conditions,  the  effects  of  axial  and  lateral  loads  may 
be  treated  independently. 


Axial  Behavior 

A  compressive  load  applied  to  the  head  (top)  of  the  pile  is  transferred  to  the 
surrounding  soil  by  a  combination  of  skin  friction  along  the  embedded  length 
and  end  bearing  at  the  tip  (bottom)  of  the  pile.  For  relatively  short  piles,  only 
the  end  bearing  effect  is  significant.  For  relatively  long  piles  in  soil  (excluding 
tip  bearing  piles  on  rock),  the  predominant  load  transfer  is  due  to  skin  friction. 
Unless  special  mechanical  provisions  are  present  (e.g.,  an  underreamed  tip), 
axial  tension  load  is  resisted  only  by  skin  friction. 
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Lateral  Behavior 


Piles  are  often  required  to  support  loads  applied  perpendicular  to  their  longi¬ 
tudinal  axes  (lateral  loads).  As  stated  previously,  lateral  load  resistance  is 
largely  independent  of  axial  effects.  However,  a  high  axial  compression  may 
interact  with  lateral  displacements  (the  beam-column  effect)  to  increase  lateral 
displacements,  bending  moments,  and  shears. 


Battered  Piles 

If  the  horizontal  loads  imparted  to  the  pile  foundation  are  large,  a  foundation 
consisting  solely  of  vertical  piles  may  not  possess  sufficient  lateral  resistance. 

In  such  circumstances,  battered  (inclined)  piles  are  installed  to  permit  the  hori¬ 
zontal  foundation  load  to  be  supported  by  a  component  of  the  axial  pile/soil 
resistance  in  addition  to  the  lateral  resistance. 


Classical  Analysis  and/or  Design  Procedures 

Single  piles 

Prior  to  the  development  of  reliable  computer  programs,  the  design  of  a 
single  pile  was  based  primarily  on  the  ultimate  load  capacity  of  the  pile  as  deter¬ 
mined  from  a  load  test  or  from  semi-empirical  equations.  The  allowable  or 
working  load  to  which  the  pile  could  be  subjected  was  taken  as  some  fraction  of 
the  ultimate.  Little,  if  any,  emphasis  was  placed  on  the  load-displacement 
behavior  of  the  pile.  Design  methodology  used  in  the  Corps  of  Engineers  is 
documented  in  Engineer  Manual  1110-2-2906  (U.S.  Army  Corps  of  Engineers 
(USACE)  1995). 

Pile  groups 

Classical  methods  (e.g.,  Culmann’s  method,  the  Common  Analytical  Method, 
the  Elastic  Center  Method,  the  Moment-of-Inertia  Method,  etc.)  of  analysis  for 
pile  groups  were  based  on  numerous  simplifying  assumptions  to  allow  the 
numerical  calculations  to  be  performed  by  hand.  Common  to  these  methods  are 
the  assumptions  that  only  the  axial  resistance  of  the  piles  is  significant  and  that 
the  pile  cap  is  rigid.  Force  and  moment  equilibrium  equations  are  used  to  allot 
the  foundation  loads  to  the  individual  piles.  No  attempt  is  made  in  these  meth¬ 
ods  to  consider  force-displacement  compatibility  (the  soil-structure  interaction 
effect).  It  has  been  shown  that  these  classical  methods  frequently  result  in 
unconservative  designs. 
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State-of-the-Corps-Art  Methods 
for  Hydraulic  Structures 

System  modelling 

Rational  designs  must  be  based  on  solutions  in  which  equilibrium  and  force- 
displacement  compatibility  are  simultaneously  satisfied.  Ongoing  research  has 
resulted  in  the  development  of  mathematical  models  for  the  pile/soil  system 
which  permit  analysis  of  the  entire  range  of  load-displacement  response  for 
single  piles  subjected  to  axial  and/or  lateral  loads.  Methods  have  been 
developed  for  the  design  of  pile  groups  in  which  the  soil-structure  interaction 
characteristics  of  single  piles  have  been  incorporated.  These  methods  and  the 
considerations  leading  to  their  development  are  described  in  detail  in  Chapters 
2-4.  A  synopsis  is  provided  in  the  following  paragraphs. 


Axially  loaded  piles 

For  analysis  of  a  pile  subjected  to  axial  loads,  the  soil  surrounding  the  em¬ 
bedded  length  of  the  pile  is  modelled  as  a  distribution  of  springs  which  resist 
longitudinal  displacements  of  the  pile.  The  resistance  of  the  soil  springs  is  rep¬ 
resentative  of  the  skin  friction  of  the  soil  on  the  pile.  The  effect  of  tip  resistance 
is  represented  by  a  concentrated  spring.  The  characteristics  of  these  springs  are 
provided  in  the  form  of  resistance-displacement  (load-transfer)  curves  represent¬ 
ing  the  skin  friction  effects  (Seed  and  Reese  (1957),  and  other  references)  and  a 
force-displacement  curve  representing  the  tip  reaction. 

The  load-transfer  curves  and  tip  reaction  curves  have  been  obtained  from 
field  tests  of  instrumented  piles  subjected  to  axial  compression.  Research  is 
continuing  to  permit  evaluation  of  load-transfer  curves  for  piles  in  tension.  The 
underlying  principles  on  which  the  load-transfer  curves  and  tip  reaction  curve 
are  based  and  the  modelling  of  the  pile/soil  system  are  presented  in  Chapter  2. 


Laterally  loaded  piles 

The  soil  which  resists  displacements  of  a  laterally  loaded  pile  is  also  replaced 
by  distributed  springs.  The  force-displacement  characteristics  of  the  springs  are 
presented  as  curves  which  have  been  extracted  from  field  tests  of  laterally 
loaded  piles.  Techniques  for  lateral  load  analysis  are  discussed  in  Chapter  3. 


Pile  head  stiffnesses 

Computer  programs  (e.g.  CAXPILE,  CPGS,  COM624G)  are  available  which 
permit  the  analysis  of  load-displacement  response  of  a  pile/soil  system  up  to  an 
ultimate  or  failure  condition.  The  relationship  between  load  and  displacement 
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tends  to  be  essentially  linear  through  the  range  of  loads  usually  allowed  (the 
working  loads)  in  design.  The  relationship  becomes  highly  nonlinear  as  an 
ultimate  condition  is  neared.  For  design  purposes,  the  linearly  elastic  relation¬ 
ship  between  head  loads  and  head  displacements  is  usually  presented  as  a  matrix 
of  stiffness  coefficients.  These  coefficients  may  be  extracted  from  the  full  range 
analyses  for  axially  or  laterally  loaded  piles  cited  above.  In  addition,  the  stiff¬ 
ness  coefficients  may  be  estimated  using  linearized  solutions.  These  processes 
are  discussed  in  Chapters  2  and  3. 


Pile  groups 

Pile  group  behavior  is  analyzed  by  the  procedure  suggested  by  Saul  (1968). 
The  method  considers  both  equilibriun  and  force-displacement  compatibility  in 
distributing  the  loads  on  the  foundation  among  the  individual  piles.  The  process 
requires  an  evaluation  of  the  linearized  pile  head  stiffness  matrix  for  each  pile  in 
the  group.  The  pile  head  stiffness  matrix  may  be  evaluated  by  the  single  pile 
analysis  procedures  alluded  to  above.  However,  the  evaluation  must  account  for 
the  effects  of  the  proximity  of  adjacent  piles. 

Although  the  group  analysis  method  was  originally  developed  for  linear  sys¬ 
tems  with  rigid  pile  caps,  it  has  been  extended  to  allow  for  flexible  caps  and,  by 
iterative  solutions,  can  account  for  nonlinear  behavior  (e.g.  CPGA).  The  method 
is  described  in  detail  in  Chapter  4. 
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2  Single  Axially  Loaded  Pile 
Analysis 


Introduction 

A  schematic  of  an  axially  loaded  pile  is  shown  in  Figure  1.  In  the  discussions 
which  follow,  the  pile  is  assumed  to  be  in  contact  with  the  surrounding  soil  over 
its  entire  length.  Consequently,  the  embedded  length  and  the  total  length  of  the 
pile  are  the  same.  The  effect  of  a  free-standing  portion  of  the  pile  will  be 
discussed  later. 

The  pile  is  assumed  to  have  a  straight  centroidal  axis  (the  z-axis,  positive 
downward)  and  is  subjected  to  a  centric  load  at  the  head  (top  of  the  pile)  P 0. 
Displacements  parallel  to  the  axis  of  the  pile  are  denoted  w  and  are  positive  in 
the  positive  z-direction.  The  pile  material  is  assumed  to  be  linearly  elastic  for  all 
levels  of  applied  loads.  “Ultimate”  conditions  referred  to  subsequently  indicate 
that  a  limit  has  been  reached  in  which  any  additional  head  load  would  cause 
excessive  displacements. 

The  major  research  efforts  devoted  toward  investigation  of  axially  loaded 
piles  have  been  performed  for  homogeneous  soil  media.  Only  in  limited  cases 
has  the  effect  of  nonhomogeneity  been  considered.  In  most  cases  the  effects  of 
layering  in  the  soil  profile  and/or  lateral  variations  in  soil  characteristics  can 
only  be  approximated. 


Load-Transfer  Mechanism 

The  head  load  P0  is  transferred  to  the  surrounding  soil  by  shear  stresses  (skin 
friction)  along  the  lateral  pile/soil  interface  and  by  end-bearing  at  the  pile  tip 
(bottom  of  the  pile).  The  rate  at  which  the  head  load  is  transferred  to  the  soil 
along  the  pile  and  the  overall  deformation  of  the  system  are  dependent  on 
numerous  factors.  Among  these  are:  (a)  the  cross  section  geometry,  material, 
length,  and,  to  a  lesser  extent,  the  surface  roughness  of  the  pile;  (b)  the  type  of 
soil  (sand  or  clay)  and  its  stress-strain  characteristics;  (c)  the  presence  or  absence 
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Figure  1.  Axially  loaded  pile 


of  groundwater;  (d)  the  method  of  installation  of  the  pile;  and,  (e)  the  presence 
or  absence  of  residual  stresses  as  a  result  of  installation. 

A  heuristic  approach  has  been  followed  to  reduce  the  complex  three-  dimen¬ 
sional  problem  to  a  quasi  one-dimensional  model  (illustrated  in  Figure  2)  which 
is  practicable  for  use  in  a  design  environment.  In  the  one-dimensional  model, 
the  soil  surrounding  the  pile  is  replaced  by  a  distribution  of  springs  along  the 
length  of  the  pile  and  by  a  concentrated  spring  at  the  pile  tip  which  resist  axial 
displacements  of  the  pile.  The  characteristics  of  these  springs  are  presented  in 
the  form  of  curves  which  provide  the  magnitude  of  unit  skin  friction  (f-w  curves) 
or  unit  tip  reaction  ( q-w  curve)  as  a  function  of  pile  displacement.  The  nomen¬ 
clature  used  to  define  axial  curves  is  based  on  unit  skin  friction  f \  unit  tip 
reaction  q,  and  w  =  displacement  in  the  z-direction  for  axial  loads. 
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Figure  2.  One-dimensional  model  of  axially  loaded  pile 


Th ef-w  and  q-w  curves  have  been  developed  using  the  principles  of  contin¬ 
uum  and  soil  mechanics  and/or  from  correlations  with  the  results  of  field  tests  on 
instrumented  axially  loaded  piles.  Several  different  criteria  are  presented  below 
for  development  of f-w  and  q-w  curves.  The  reliability  of  any  method  in 
predicting  the  behavior  of  a  particular  pile  depends  on  the  similarity  of  the 
system  under  investigation  with  the  database  used  to  establish  the  method.  Most 
of  the  methods  account  explicitly  or  implicitly  for  the  three  factors  cited  on 
page  5  (a,  b,  and  c).  In  all  cases  the  pile  is  assumed  to  be  driven  into  the  soil  or 
to  be  a  cast-in-place  pier.  Only  one  of  the  procedures  attempts  to  account  for  the 
effects  of  residual  stresses;  the  remaining  methods  exclude  these  effects. 
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Synthesis  of  f-w  Curves  for  Piles  in 
Sand  Under  Compressive  Loading 

Mosher  (1984) 

Mosher  (1984)  utilized  the  results  of  load  tests  of  prismatic  pipe  piles  driven 
in  sand  and  the  work  of  Coyle  and  Castello  (1981)  to  arrive  at  the  hyperbolic 
representation  of  the  f-w  curve  (see  Figure  3). 


Table  1 

kf  (psf/in.)  as  Function  of  Angle  of  Internal  Friction  of  Sand  for 
Method  SSF1 

Angle  of  Internal  Friction  (degrees)1 

kf  (psf/in.) 

28-31 

6,000-10,000 

32-34 

10,000-14,000 

35-38 

14,000-18,000 

1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units  is  presented 

||  on  page  x. 

Unit  Ultimate  Side  Friction  -  f  (tsf) 
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Figure  4.  Ultimate  side  friction  for  Method  SSF1 
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The  effects  of  groundwater,  layering,  and  variable  pile  diameter  may  be 
accounted  for  approximately  by  adjusting  the  relative  depth  at  each  point  as  fol¬ 
lows.  The  effective  depth  z'  below  ground  surface  is  obtained  by  dividing  the 
effective  vertical  soil  pressure  at  a  point  by  the  effective  unit  weight  at  that 
point;  the  relative  depth  is  obtained  by  dividing  the  effective  depth  by  the  pile 
diameter  at  that  point.  This  approximation  will  result  in  unrealistic  discon¬ 
tinuities  in  the  distribution  of  f-w  curves  at  soil  layer  boundaries,  at  the  location 
of  a  subsurface  groundwater  level,  and  at  changes  in  pile  diameter. 

The  method  may  also  be  extended  to  approximate  the  behavior  of  noncircular 
cross  sections  using  the  equivalent  radius  of  the  pile  as  indicated  in  Figure  5. 


Kraft,  Ray,  and  Kagawa  (1981) 

Numerous  analyses  (Randolph  and  Wroth  1978;  Vesic  1977;  Kraft,  Ray,  and 
Kagawa  1981;  Poulos  and  Davis  1980)  have  been  performed  in  which  the 
pile/soil  system  is  assumed  to  be  radially  symmetric  and  the  soil  is  assumed  to 
be  a  vertically  and  radially  homogeneous,  elastic  medium.  The  principles  of 
continuum  mechanics  as  well  as  finite  element  methods  have  been  used  to  arrive 
at  the  relationship  between  side  friction  and  axial  pile  displacement.  The 
process  due  to  Kraft,  Ray,  and  Kagawa  (1981)  is  outlined  below. 

Shear  stresses  are  assumed  to  decay  radially  in  the  soil  according  to 

t  =  fRJr  (2) 
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where 


t  =  shear  stress  in  the  soil 
/  =  side  friction  at  pile/soil  interface 
R  =  pile  radius 

r  =  radial  distance  from  the  pile  centerline 

If  radial  deformations  of  the  soil  are  ignored,  the  shear  strain  at  any  point  in 
the  soil  may  be  expressed  as 

Y  =  *!L  =  1-  =  (3) 

dr  G  Gr 


where  G  is  the  soil  shear  modulus  of  elasticity. 

The  axial  displacement  at  the  interface  is  obtained  by  integrating  Equation  3 
to  obtain 


w 


=  fR 

G 


r}*L-„fAin 

\ 

rm 

J  r  G 

where 


(4) 


w  =  axial  displacement  of  the  pile 

rm  =  a  limiting  radial  distance  beyond  which  deformations  of  the  soil  mass  are 
negligible 

Randolph  and  Wroth  suggested 

rm  =  2Zp(l  -  v)  (5) 

where 


L  -  embedded  length  of  the  pile 

p  =  a  factor  to  account  for  vertical  nonhomogeneity  of  the  soil  medium  to  be 
discussed  later 


v  =  Poisson's  ratio  for  the  soil 

Combination  of  Equations  4  and  5  yields  a  linear  relationship  between  pile 
displacements  and  side  friction  as 
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fR  , 

w  -  - —  In 
G 


(6) 


2Ip(l  -  v) 

R 

The  side  friction-displacement  relationship  expressed  in  Equation  6  is 
appropriate  only  for  very  small  displacements.  To  account  for  deviations  from 
linearity,  a  hyperbolic  variation  in  side  friction-displacements  proposed  by 
Kraft,  Ray,  and  Kagawa  (1981)  is 


2Lp(l  -  v) 
R 


B 

f 

J max 


l 


/ 

J  TTU 


(?) 


where  Rf  is  a  curve  fitting  parameter  (Kraft,  Ray,  and  Kagawa  1981)  which  may 
be  taken  as  0.9  for  most  conditions.  The  value  of fmax  may  be  obtained  from  the 
curves  due  to  Mosher  (Figure  4)  or  may  be  estimated  as  suggested  under  method 
SSF3  which  follows. 

After  fmax  has  been  reached,  the  f-w  curve  becomes  a  horizontal  line  at  fmax  . 

Th ef-w  curve  produced  by  this  method  is  illustrated  in  Figure  6  by  the  solid 
curve  0  -  1  -  2  . 

Some  soils  exhibit  a  degradation  in  strength  after  a  maximum  resistance  has 
been  reached.  The  results  of  a  direct  shear  test  for  a  softening  soil  illustrated  in 
Figure  7  are  used  to  construct  the  descending  branch  of  the  f-w  curve  for 
softening  soils  shown  by  the  dashed  line  in  Figure  6  as  follows.  The  displace¬ 
ment  beyond  the  maximum  fmax  required  to  reduce  the  side  friction  to  its  residual 
value  is  obtained  by  adjusting  the  direct  shear  displacement  to  account  for  elastic 
rebound  of  the  pile  due  to  the  reduction  in  side  friction.  This  adjustment  is  given 
by 


w 


frnJj 

G 


2Lp(l  -  v)  R 
R  f 

2Zp(l  —  V)  f resicftf 

In 

-  In 

R  f 

J  max 

l-Rf 

j  f resicft'f 

f 

J  max 

(8) 


The  softening  portion  of  the  f-w  curve  is  obtained  by  scaling  the  normalized 
direct  shear  curve  to  the  f-w  curve  (dashed  line  1-3  in  Figure  6). 

Because  the  shear  strength  of  sands  increases  with  depth  (i.e.  confining 
pressure),  the  shear  modulus  G  is  not  constant  along  the  length  of  the  pile. 

Finite  element  analyses  have  indicated,  for  a  linear  increase  in  G  with  depth,  the 
value  of  p  to  be 
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Displacement  —  w 


Figure  6.  f-w  curve  by  Method  SSF2 


P 


(9) 


where 

Gm  =  soil  shear  modulus  at  mid-depth  of  the  pile 
G,  =  shear  modulus  at  the  pile  tip 

The  preceding  equations  also  assume  that  the  soil  modulus  G  is  unaffected 
by  the  pile  installation.  Randolph  and  Wroth  (1978)  performed  finite  element 
analyses  for  two  hypothetical  variations  of  shear  modulus  radially  away  from  the 
pile.  These  variations  and  the  effective  shear  modulus  were: 

a.  G  =  GV4  for  1  <  r/R  <  1.25  ;  G  =  for  r/R  >  1.25  which  produced 


G°° 


3  In  1.25  +  In 


2  £p(l  -  v) 
R 


(10) 
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Figure  7.  Direct  shear  test  of  softening  soil 

b.  G  =  GW4  for  1  <  r/R  <  1.25  ;  G  =  for  r/R.2  ;  G  varied  linearly 
between  1.25  <  r/R  <2  which  produced 


Geff  = 


Gx 


5  In  1.25  +  In 


2Z,p(l  -  v) 


R 


where 


(11) 


Geff=  reduced  effective  shear  modulus 
=  shear  modulus  of  the  undisturbed  soil 

The  shear  modulus  G  in  the  preceding  equations  should  be  evaluated  at  a  low 
strain  value  such  as  in  the  range  of  values  obtained  from  seismic  velocity  tests 
conducted  in  situ  or  from  resonant  column  tests.  As  an  alternative,  the  following 
expression  may  be  used 

G  =  83k  v/o^  (12) 
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where 


k=  a  function  of  relative  density,  varying  from  50  at  a  relative  density  of 
60  percent  to  70  at  a  relative  density  of  90  percent 

o'0  =  mean  effective  stress  in  the  soil  (vertical  stress  plus  two  times  horizontal 
stress);  with  G  and  a'0  in  psi 

Vijayvergiya  (1977) 

Vijayvergiya  (1977)  proposed  a  relationship  between  side  friction  and  pile 
displacement  of  the  form 


where  wc  is  the  displacement  required  to  develop/^.  For  w  greater  than  wc  ,/ 
remains  constant  at  fmax.  Vijayvergiya  gives  limiting  values  of  fmax  as  1  tsf  for 
clean  medium  dense  sand,  0.85  tsf  for  silty  sand,  0.7  tsf  for  sandy  silt,  and  0.5  tsf 
for  silts.  The  suggested  values  of  wc  range  from  0.2  to  0.3  in.  for  nominal  sized 
piles.  A  typical  f-w  curve  by  this  method  is  shown  as  the  solid  curve  in  Figure  8. 


Figure  8.  f-w  curve  by  Method  SSF3 
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Coyle  and  Sulaiman  (1967) 

Coyle  and  Sulaiman  (1967)  performed  tests  on  miniature  piles  in  sand  and 
correlated  the  laboratory  results  with  data  from  field  tests  of  instrumented  piles 
in  sand.  They  concluded  that  skin  friction  increases  with  pile  deflection  up  to 
pile  displacements  of  0.1  to  0.2  in.  They  further  concluded  that  the  ratio  of  skin 
friction  to  soil  shear  strength  is  high  (greater  than  one)  near  the  ground  surface 
and  decreases  to  a  limiting  value  of  0.5  with  increasing  depth.  Two  curves,  as 
shown  in  Figure  9,  were  proposed  for  the  analysis  of  axially  loaded  piles  in  sand. 
Curve  A  was  proposed  for  use  at  depths  less  than  20  ft  below  the  surface  and 
Curve  B  for  depths  greater  than  20  ft. 


Figure  9.  f-w  curves  by  Method  SSF4 


Briaud  and  Tucker  (1984) 

Analyses  using  the/-w  curves  discussed  above  do  not  consider  the  presence 
of  residual  stresses  in  the  pile/soil  system  which  result  from  the  installation 
process.  Field  tests  of  instrumented  piles  indicate  that  significant  residual 
stresses  may  be  encountered  in  long,  flexible  piles  driven  in  sand  or  gravel  (see, 
for  instance,  Mosher  (1984)).  If  the  f-w  curves  and  tip  reaction  representation 
(see  later)  are  both  based  on  ignoring  residual  stresses,  the  predicted  pile  head 
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displacement  at  any  load  will  be  essentially  unaffected.  However,  the  distribu¬ 
tion  of  axial  load  and  the  predicted  tip  reaction  may  be  in  error.  Briaud  and 
Tucker  (1984)  extracted  the  residual  stresses  from  field  tests  of  piles  in  sand.  A 
hyperbolic  representation  of  the  f-w  curve  (Figure  10)  was  proposed  for  inclu¬ 
sion  of  the  effects  of  residual  stresses  as 


/=- - - -  -/v  ' 

1  w 

+ 

(14) 

kf  fmax  fr 

with 

kf  =  5.01  (/V/27 

(15) 

fmax  -  0.224fVjr 

(16) 

KD 

n 

(17) 

fr  =  5.57IQ  4  <fmax 

A* 

(18) 

where 


Ns  =  number  of  blows  per  foot  in  a  standard  penetration  test 
C  =  21 IR  =  pile  circumference 
A  =  pile  cross  section  area 
E  =  pile  modulus  of  elasticity 
L  =  length  of  pile 
At  =  tip  reaction  area 

As  =  CL  =  area  of  pile-soil  interface;  with  kf  in  tsf/in.;  and  fmax  and/r  in  tsf 

Synthesis  of  f-w  Curves  for  Piles  in 
Clay  Under  Compressive  Loading 

Coyle  and  Reese  (1 966) 

The  results  of  load  tests  of  instrumented  piles  in  clay  as  well  as  the  results  of 
laboratory  tests  of  model  pile/soil  systems  were  used  by  Coyle  and  Reese  (1966) 
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Figure  10.  f-w curve  by  Method  SSF5 


to  establish  the  three  load  transfer  curves  shown  in  Figure  1 1 .  Curve  A  is 
applicable  for  points  along  the  pile  from  the  ground  surface  to  a  depth  of  10  ft, 
curve  B  applies  for  depths  from  10  ft  to  20  ft,  and  curve  C  is  applicable  for  all 
depths  below  20  ft. 

The  relationship  between  maximum  side  friction  and  soil  shear  strength 
provided  by  Coyle  and  Reese  is  shown  in  Figure  12. 


Aschenbrener  and  Olson  (1984) 

Data  obtained  from  a  large  number  of  field  load  tests  of  piles  in  clay  were 
examined  by  Aschenbrener  and  Olson  (1984)  with  the  intent  to  devise  load 
transfer  relationships  which  provided  the  best  fit  to  the  diverse  pile  and  soil 
properties  represented  by  the  database.  The  simple  bilinear  relationship  shown 
in  Figure  13  was  selected  as  a  result  of  their  study. 

Aschenbrener  and  Olson  expressed  the  relationship  between  fmax  and  soil 
shear  strength  as 


(19) 
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1.0 


Curve  B 


Curve  A 


Displacement 


Figure  1 1 .  f-w curves  by  Method  CSF1 


Figure  12.  Side  friction  -  soil  strength  relation  for  Method  CSF1 
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Figure  13.  f-w curve  by  Method  CSF2 
where 

a  =  a  proportionality  factor 
su  =  undrained  shear  strength 

Aschenbrener  and  Olson  were  able  to  evaluate  a  from  the  field  test  data  as 

a  _  (20) 

SuAs 

where 


Pou  =  pile  head  load  at  failure 
Pm  =  tip  load  at  failure 
su  =  undrained  shear  strength 
As  =  area  of  pile-soil  interface 

In  a  design  situation,  the  ultimate  head  and  tip  loads  will  not  be  known.  For 
design,  the  value  of  a  may  be  obtained  from  the  curves  provided  by  Semple  and 
Rigden  (1984)  shown  in  Figure  14  as 

a  =  ajj,  (21) 
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where 


ap  =  peak  strength  reduction  factor  from  Figure  14a 
a,  =  length  factor  from  Figure  14b 


In  Figure  14,  su  is  the  undrained  shear  strength;  ov  is  the  effective  overburden 
pressure;  L  is  the  length  of  pile;  and,  R  is  the  pile  radius. 


Kraft,  Ray,  and  Kagawa  (1981) 

The  procedure  of  Method  SSF2  due  to  Kraft,  Ray,  and  Kagawa  (1981) 
described  previously  for  sand  side  friction,  may  be  applied  to  piles  in  clay.  For 
clays,  the  shear  modulus  may  again  be  evaluated  from  seismic  tests,  from 
resonant  column  tests,  approximated  as  400  to  500  times  su ,  or  evaluated  from 
the  modulus  of  elasticity  as  El 3  for  undrained  conditions  and  EH. 15  for  drained 
conditions. 


Chapter  2  Single  Axially  Loaded  Pile  Analysis 


21 


Heydinger  and  O’Neill  (1986) 


Finite  element  and  finite  difference  analyses  were  performed  by  Heydinger 
and  O'Neill  (1986)  to  develop  f-w  curves  for  piles  in  clay.  An  axisymmetric 
model  including  interface  elements  to  account  for  slippage  of  the  pile-soil 
system  was  used  in  the  finite  element  analyses.  An  unconsolidated-undrained 
condition  was  assumed  to  exist  in  the  .soil  and  the  initial  mean  effective  stresses 
were  computed  from  radial  consolidation  theory  in  which  the  pile  installation 
process  was  represented  by  an  expanding  cylindrical  cavity.  A  general  equation 
for  the  f-w  curves  (illustrated  in  Figure  15)  was  selected  as 


w 
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bf 

X  w  i 


II  lm 


(22) 


where  the  parameters  E{  and  m  were  determined  by  statistical  correlations  with 
the  analytical  results  as 


Ef  =  jE/exp 


0.36  +  0.38  In  — 
2R  J 


(23) 


and 


m  -  exp 


0.12  +  0.54  In  ^  -  0.42  In  — 
Pa  2R 


(24) 


where 

E  =  initial  undrained  modulus  of  elasticity  of  the  soil  at  the  depth  of  interest 

Eavg  =  the  average  initial  undrained  modulus  of  elasticity  over  the  entire  length 
of  the  pile 

pa  =  atmospheric  pressure  in  the  same  units  as  Eavg 

The  value  off  should  be  measured  at  very  low  strains.  An  approximation  for  E 
is  cited  as  1,200  to  1,500  times  su . 


Vijayvergiya  (1977) 

Vijayvergiya  (1977)  indicated  that  Equation  13  for  f-w  curves  in  sand, 
Method  SSF3,  is  applicable  for  piles  in  clay.  As  for  piles  in  sand,  Vijayvergiya 
suggests  values  of  the  critical  pile  displacement  wc  of  0.2  to  0.3  in.  Although  a 
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Displacement  Ratio  -  w/2R 


Figure  15.  f-w  curve  by  Method  CSF4 

method  for  evaluating/^  is  presented  by  Vijayvergiya,  he  suggests  that  other  less 
complex  methods  are  equally  suitable,  e.g.,  the  process  of  Method  CSF2 
discussed  previously. 


Tip  Reactions 

The  influence  of  the  tip  reaction  on  the  axial  load-displacement  behavior 
depends  on  the  relative  stiffness  of  the  pile  as  well  as  side  friction  stiffness  of 
the  soil.  In  the  following  paragraphs  several  curves  are  presented  for  assessing 
the  tip  reaction  as  a  function  of  the  tip  displacement.  In  general  these  curves 
have  been  developed  primarily  from  a  consideration  of  the  properties  of  the  soil 
at  the  tip  elevation.  However,  numerous  theoretical  studies  (see,  for  instance, 
Randolph  and  Wroth  (1978))  have  indicated  that  the  tip  reaction  depends  on  the 
characteristics  of  the  soil  both  above  and  below  the  tip  elevation.  Some  of  the 
methods  for  developing  q-w  curves  for  the  tip  reaction  account  for  the  profile  in 
the  vicinity  of  the  tip  by  using  average  soil  properties.  Other  methods,  derived 
from  test  results  where  the  soil  at  the  test  site  was  relatively  homogeneous,  are 
dependent  on  the  properties  of  the  soil  at  the  tip. 

The  curves  presented  below  are  for  unit  tip  reaction  (i.e.  force  per  unit  of  tip 
area).  To  evaluate  the  total  tip  reaction,  this  unit  force  must  be  multiplied  by  the 
area  of  the  pile  tip  actually  bearing  on  the  soil.  For  solid  or  closed-end  piles  the 
tip  bearing  area  is  reasonably  taken  as  being  equal  to  the  gross  cross  section 
area.  For  open-end  piles  (e.g.  pipes)  or  H-piles  the  effective  tip  area  may  be  as 
little  as  the  material  area  of  the  pile  or  may  be  as  much  as  the  gross  section  area. 
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Synthesis  of  q-w  Curves  for  Piles  in 
Sand  Under  Compressive  Loading 

Mosher  (1984) 

Mosher  (1984)  expanded  the  work  of  Coyle  and  Castello  (1981)  to  determine 
the  q-w  relationship  for  piles  in  sand.  Mosher  proposed  the  exponential  q-w 
curve  shown  in  Figure  16.  Values  of  ultimate  unit  tip  reaction  qmax  are  given  as  a 
function  of  relative  depth  (L/2R)  in  Figure  17. 


Kraft,  Ray,  and  Kagawa  (1981) 

Kraft,  Ray,  and  Kagawa  (1981)  did  not  attempt  to  produce  a  q-w  curve 
corresponding  to  their  analytical  f-w  curve,  but  approximated  the  q-w  relation¬ 
ship  by  the  elastic  solution  for  a  rigid  punch  according  to 

w  =  t1  -^1  I  (25) 

E  ' 

where 

w  =  tip  displacement 
R  =  radius  of  pile  tip  reaction  area 
q  =  tip  pressure 

v  =  Poisson's  ratio  for  the  soil  at  the  tip 

E  =  secant  modulus  of  elasticity  of  the  soil  appropriate  to  the  level  of  soil 
stress  associated  with  q 

I,  -  influence  coefficient  ranging  from  0.5  for  long  piles  to  0.78  for  very  short 
piles 


Vijayvergiya  (1977) 

Vijayvergiya  (1977)  proposed  an  exponential  representation  for  the  q-w 
curve  for  a  pile  in  sand  similar  to  those  in  Method  ST1.  For  w  <wc 


24 


Chapter  2  Single  Axially  Loaded  Pile  Analysis 


Figure  1 6.  q-w  curve  by  Method  ST  1 
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Figure  17.  Ultimate  tip  resistance  for  Method  SF1 
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where  wc  is  the  critical  tip  displacement  given  by  Vijayvergiya  as  ranging  from 
3  to  9  percent  of  the  diameter  of  the  tip  reaction  area.  For  w>  wc,q  =  qmax. 
Vijayvergiya  did  not  suggest  adjusting  the  exponent  to  account  for  density. 


Briaud  and  Tucker  (1984) 

Briaud  and  Tucker  (1984)  offer  a  means  of  accounting  for  the  presence  of 
residual  stresses  due  to  pile  installation  on  the  tip  reaction.  The  hyperbolic 
relationship  between  unit  tip  reaction  and  tip  displacement  shown  in  Figure  18  is 
given  by 


II 

h 

+ 

(27) 

kq  Qmax  Qr 

-  19-75(1 V)“ 

(28) 

qr  =  5.57IQ 

(29) 

Kg  =  467.1  (TV)00065 

(30) 

where 

N  =  uncorrected  average  blow  count  of  a  standard  penetration  test  over  a 
distance  of  four  diameters  on  either  side  of  the  tip 

kq  =  initial  slope  of  the  q-w  curve  in  tsf/in. 

Qmax  >  Qr =  ultimate  and  residual  unit  tip  resistances,  respectively,  in  tsf.  Other 
terms  are  defined  on  page  6 


Coyle  and  Castello  (1981) 

Coyle  and  Castello  (1981)  provided  ultimate  tip  reactions  based  on  correla¬ 
tions  for  instrumented  piles  in  sand  as  shown  in  Figure  19.  Coyle1  recom¬ 
mended  the  tip  reaction  curve  shown  in  Figure  20. 


1  Unpublished  Class  Notes,  1977,  H.  M.  Coyle,  “Marine  Foundation  Engineering,”  Texas  A&M 
University,  College  Station,  TX. 
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Figure  18.  q-w  curve  by  Method  SF4 

Synthesis  of  q-w  Curves  for  Piles  in 
Clay  Under  Compressive  Loading 

Aschenbrener  and  Olson  (1984) 

Data  for  tip  load  and  tip  settlement  were  not  recorded  in  sufficient  detail  in 
the  database  considered  by  Aschenbrener  and  Olson  (1984)  to  allow  establishing 
a  nonlinear  q-w  relationship.  It  was  concluded  that  the  sparsity  and  scatter  of 
field  data  warranted  nothing  more  complex  than  a  simple  elasto-plastic  relation¬ 
ship.  In  their  representation,  q  varies  linearly  with  w  reaching  qmax  at  a  displace¬ 
ment  equal  to  1  percent  of  the  tip  diameter  and  remains  constant  at  qmax  for  larger 
displacements.  Ultimate  tip  reaction  was  evaluated  according  to 


where 

su  =  undrained  shear  strength 
Nc  =  bearing  capacity  factor 
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Figure  1 9.  Ultimate  tip  resistance  for  Method  ST5 


Test  data  indicated  that  Nc  varied  from  0  to  20  and  had  little  correlation  with 
shear  strength.  When  ultimate  tip  reaction  was  not  available  from  recorded  data, 
Aschenbrener  and  Olson  used  a  conventional  value  for  Nc  equal  to  9. 


Vijayvergiya  (1977) 

Vijayvergiya  (1977)  recommends  that  the  exponential  q-w  curve  for  sand  as 
discussed  on  pages  24-26  is  applicable  for  piles  in  clay.  He  indicates  that  qmax 
can  be  calculated  from  Equation  3 1  above  but  provides  no  guidance  for  the 
selection  of  Nc. 


Other  Considerations 

Uplift  loading 

For  some  design  cases  it  may  be  necessary  to  evaluate  the  behavior  of  an 
axially  loaded  pile  for  uplift  (tension)  loading.  Considerably  less  is  known  about 
uplift  loading  than  about  compression  loading.  However,  it  is  believed  to  be 
sufficiently  accurate  to  analyze  prismatic  piles  in  clay  under  uplift  using  the  same 
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Axial  Displacement  -  w  (in.) 


Figure  20.  q-w  curve  by  Method  SF5 


procedures  used  for  compression  loading,  except  that  the  tip  reaction  should  be 
omitted  unless  it  is  explicitly  accounted  for  as  discussed  below.  In  sands,  use  of 
the  same  procedures  employed  in  compression  loading  is  recommended,  with 
the  exception  that  should  be  reduced  to  70  percent  of  the  maximum  compres¬ 
sion  value. 


For  the  methods  that  explicitly  include  residual  driving  stress  effects  in 
nonlinear/w  and  q-w  curves  (pages  16-17  and  26),  it  is  recommended  that  the 
appropriate  curves  for  uplift  loading  be  generated  by  extending  the  solid  curves 
in  Figures  10  and  18  in  the  negative  w  direction  with  the  same  initial  slopes  as 
exist  in  the  positive  w  direction  and  assuring  that  the  q-w  curve  terminates  at  q  = 
0.  That  is 

/=- - - - 

I _ — _  (32; 

kf  Q-Vmax  +  fr 

where  w  is  negative  andfmax,fr,  and  kf  are  positive.  And 


9  = 


9r 


(33) 
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where  w  is  negative  and  qr  and  kq  are  positive.  All  parameters  appearing  in 
Equations  32  and  33  are  evaluated  as  for  compressive  loading. 


Bearing  on  Rock 

The  tip  reaction-tip  displacement  relationship  for  a  pile  driven  to  bearing  on 
rock  may  be  assumed  to  be  linear.  The  tip  reaction  stiffness  given  by  Equation 
25  may  be  used  where  the  modulus  of  elasticity  and  Poisson's  ratio  should  reflect 
the  characteristics  of  the  surficial  zone  of  the  rock.  The  influence  coefficient  /, 
in  Equation  25  may  be  taken  as  0.78  for  very  sound  rock  but  should  be  reduced 
to  account  for  such  effects  as  fracturing  of  the  rock  surface  due  to  driving. 


Cyclic  Loading 

Studies  have  shown  (Poulos  1983)  that  the  principal  concern  associated  with 
cyclic  axial  loading  is  the  tendency  for  fmax  to  reduce  as  the  ratio  of  the  cyclic 
component  of  axial  load  Poc  to  the  ultimate  static  capacity  Pous  increases  beyond 
some  critical  value.  As  long  as  the  ratio  remains  below  the  failure  envelope 
shown  in  Figure  21,  no  significant  degradation  of  the  pile  capacity  or  force- 
displacement  behavior  is  likely  to  occur. 


Figure  21 .  Assessment  of  degradation  due  to  static  loading 
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Algorithm  for  Analysis  of  Axially  Loaded  Piles 

The  derivation  of  the  f-w  and  q-w  curves  from  theoretical  considerations  or 
from  experimental  data  described  in  the  preceding  sections  was  in  all  cases 
based  on  the  assumption  that  the  side  friction  f  or  tip  reaction  q  at  any  point  is  a 
function  only  of  the  pile  displacement  at  that  point  (i.e.  the  well  known  Winkler 
assumption).  For  this  assumption  and  the  one-dimensional  model  of  the  pile-soil 
system  shown  in  Figure  2,  the  governing  differential  equation  for  a  prismatic, 
linearly  elastic  pile  is 

EA  —  -  2t iRf(z,w)  =  0  (34) 

dz2 

where 

E  -  modulus  of  elasticity  of  the  pile  material 

A  =  pile  material  cross  section  area 

w  =  axial  displacement 

R  —  effective  radius  of  pile  soil  interface;  and  f[z,w)  is  the  unit  side  friction, 
which  is  a  function  of  both  position  on  the  pile  as  well  as  pile 
displacement 

Because  the  displacements  must  be  known  before  the  side  friction  fz,w)  can 
be  determined,  numerical  iterative  solutions  of  Equation  34  are  required.  The 
most  common  approach  to  the  solution  is  to  replace  the  continuous  pile-soil 
system  with  a  discretized  model  (Coyle  and  Reese  1966,  Dawkins  1982, 
Dawkins  1984)  defined  by  a  finite  number  of  nodes  along  the  pile  at  which 
displacements  and  forces  are  evaluated.  The  solution  proceeds  by  a  succession 
of  trial  and  correction  solutions  until  compatibility  of  forces  and  displacements 
is  attained  at  every  node. 


Observations  of  System  Behavior 

An  expedient  device  in  obtaining  the  numerical  solutions  described  above  is 
to  replace  the  nonlinear  f-w  and  q-w  curves  by  equivalent  linearly  elastic  springs 
during  each  iteration.  The  stiffnesses  of  these  linear  springs  are  evaluated  as  the 
secant  to  Has  f-w  or  q-w  curve  for  the  displacement  calculated  during  the 
preceding  iteration.  It  is  to  be  noticed  that  ultimate  side  friction  increases  with 
depth  while  pile  displacements  decrease  with  depth.  Hence  it  can  be  concluded 
that  the  stiffness  of  the  load  transfer  mechanism  for  side  friction  increases  with 
depth.  If  the  distribution  of  the  side  friction  for  any  given  head  load  can  be 
determined  then  a  solution  may  be  obtained  from  a  linearly  elastic  solution 
without  the  need  for  iterations. 
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3  Single  Laterally  Loaded 
Pile  Analysis 


Introduction 

Although  the  usual  application  of  a  pile  foundation  results  primarily  in  axial 
loading,  there  exist  numerous  situations  in  which  components  of  load  at  the  pile 
head  produce  significant  lateral  displacements  as  well  as  bending  moments  and 
shears.  Unlike  axial  loads,  which  only  produce  displacements  parallel  to  the 
axis  of  the  pile  (a  one-dimensional  system),  lateral  loads  may  produce  displace¬ 
ments  in  any  direction.  Unless  the  pile  cross  section  is  circular,  the  laterally 
loaded  pile/soil  system  represents  a  three-dimensional  problem.  Most  of  the 
research  on  the  behavior  of  laterally  loaded  piles  has  been  performed  on  piles  of 
circular  cross  section  in  order  to  reduce  the  three-dimensional  problem  to  two 
dimensions.  Little  work  has  been  done  to  investigate  the  behavior  of  noncircular 
cross  section  piles  under  generalized  loading.  In  many  applications,  battering  of 
the  piles  in  the  foundation  produces  combined  axial  and  lateral  loads.  However, 
the  majority  of  the  research  on  lateral  load  behavior  has  been  restricted  to 
vertical  piles  subjected  to  loads  which  produce  displacements  perpendicular  to 
the  axis  of  the  pile.  In  the  discussions  which  follow,  it  is  assumed  that  the  pile 
has  a  straight  centroidal  vertical  axis.  If  the  pile  is  nonprismatic  and  has  a 
noncircular  cross  section,  it  is  assumed  that  the  principal  axes  of  all  cross 
sections  along  the  pile  fall  in  two  mutually  perpendicular  planes  and  that  the 
loads  applied  to  the  pile  produce  displacements  in  only  one  of  the  principal 
planes. 

A  schematic  of  a  laterally  loaded  pile  is  shown  in  Figure  22.  The  x-z  plane  is 
assumed  to  be  a  principal  plane  of  the  pile  cross  section.  Due  to  the  applied 
head  shear  V0  and  head  moment  M0 ,  each  point  on  the  pile  undergoes  a  transla¬ 
tion  u  in  the  x-direction  and  a  rotation  0  about  the  y-axis.  Displacements  and 
forces  are  positive  if  their  senses  are  in  a  positive  coordinate  direction.  The 
surrounding  soil  develops  pressures,  denoted  p  in  Figure  22,  which  resist  the 
lateral  displacements  of  the  pile. 

The  principles  of  continuum  mechanics  and  correlations  with  the  results  of 
tests  of  instrumented  laterally  loaded  piles  have  been  used  to  relate  the  soil 
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Figure  22.  Laterally  loaded  pile 

lateral  resistance  p  at  each  point  on  the  pile  to  the  lateral  displacement  u  at  that 
point  (i.e.  the  Winkler  assumption).  The  relationship  between  soil  resistance  and 

lateral  displacement  is  presented  as  a  nonlinear  curve  -  the  p-u  curve.  Several 
methods  are  summarized  in  the  following  paragraphs  for  development  of  p-u 
curves  for  laterally  loaded  piles  in  both  sands  and  clays.  In  all  of  the  methods, 
the  primary  p-u  curve  is  developed  for  monotonically  increasing  static  loads. 

The  static  curve  is  then  altered  to  account  for  the  degradation  effects  produced 
by  cyclic  loads  such  as  might  be  produced  by  ocean  waves  on  offshore  struc¬ 
tures.  Methods  designated  SLAT1  and  CLAT1  through  CLAT4  have  been 
incorporated  into  the  CASE  Project  Computer  program  CPGS. 
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Load  Transfer  Mechanism  for 
Laterally  Loaded  Piles 


The  load  transfer  mechanism  for  laterally  loaded  piles  is  much  more  complex 
than  that  for  axially  loaded  piles.  In  an  axially  loaded  pile  the  axial  displace¬ 
ments  and  side  friction  resistances  are  unidirectional  (i.e.,  a  compressive  axial 
head  load  produces  downward  displacements  and  upward  side  friction  resistance 
at  all  points  along  the  pile).  Similarly,  the  ultimate  side  friction  at  the  pile-soil 
interface  depends  primarily  on  the  soil  shear  strength  at  each  point  along  the 
pile.  Because  the  laterally  loaded  pile  is  at  least  two-dimensional,  the  ultimate 
lateral  resistance  of  the  soil  is  dependent  not  only  on  the  soil  shear  strength  but 
on  a  geometric  failure  mechanism.  At  points  near  the  ground  surface  an  ultimate 
condition  is  produced  by  a  wedge  type  failure,  while  at  lower  positions  failure  is 
associated  with  plastic  flow  of  the  soil  around  the  pile  as  displacements  increase. 
In  each  of  the  methods  described  below,  two  alternative  evaluations  are  made  for 
the  ultimate  lateral  resistances  at  each  point  on  the  pile,  for  wedge  type  failure 
and  for  plastic  flow  failure,  and  the  smaller  of  the  two  is  taken  as  the  ultimate 
resistance. 


Synthesis  of  p-u  Curves 
for  Piles  in  Sand 

Reese,  Cox,  and  Koop  (1974) 

A  series  of  static  and  cyclic  lateral  load  tests  were  performed  on  pipe  piles 
driven  in  submerged  sands  (Cox,  Reese,  and  Grubbs  1974;  Reese,  Cox,  and 
Koop  1974;  Reese  and  Sullivan  1980).  Although  the  tests  were  conducted  in 
submerged  sands,  Reese  et  al.  (1980)  have  provided  adjustments  by  which  the 
p-u  curve  can  be  developed  for  either  submerged  sand  or  sand  above  the  water 
table.  The  p-u  curve  for  a  point  a  distance  z  below  the  pile  head  extracted  from 
the  experimental  results  is  shown  schematically  in  Figure  23.  The  curve  consists 
of  a  linear  segment  from  0  to  a  ,  an  exponential  variation  of  p  with  u  from  a  to 
b,  a  second  linear  range  from  b  to  c,  and  a  constant  resistance  for  displacements 
beyond  c . 

Steps  for  constructing  the  p-u  curve  at  a  depth  z  below  the  ground  surface  are 
as  follows: 

a.  Determine  the  slope  of  the  initial  linear  portion  of  the  curve  from 


kp  =  kz  (35) 

where  k  is  obtained  from  Table  2  for  either  submerged  sand  or  sand  above  the 
water  table. 
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C3  =  K  tan  <J>  tan4  P  +  tan2  (45  -  <J>/2)  (tan8  P  -  l) 


(40) 


where 

y  =  effective  unit  weight  of  the  sand 
z  =  depth  below  ground  surface  . 

K  =  horizontal  earth  pressure  coefficient  chosen  as  0.4  to  reflect  the  fact 
that  the  surfaces  of  the  assumed  failure  model  are  not  planar 

<j)  =  angle  of  internal  friction 

p  =  45  +  <{>/2 

b  =  width  of  the  pile  perpendicular  to  the  direction  of  loading 


Values  of  Cl,  C2,  C3,  and  the  depth  zcr  at  which  the  transition  from  wedge 
failure  (Equation  36)  to  flow  failure  (Equation  37)  occurs  are  shown  in 
Figure  24. 
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Figure  24.  Factors  for  calculation  of  ultimate  soil  resistance  for  laterally  loaded 
pile  in  sand 
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c.  Compute  the  lateral  resistance  for  the  transition  points  c  and  b  on  the 
curve  (Figure  23)  from 

pc  =  A  ps  (41) 

pb  =  Bps  (42> 

where  A  and  B  are  reduction  coefficients  from  Figures  25  and  26, 
respectively,  for  the  appropriate  static  or  cyclic  loading  condition.  The 
second  straight  line  segment  of  the  curve,  from  b  to  c ,  is  established  by  the 
resistances  pb  and  pc  and  the  prescribed  displacements  of  u  =  b/ 60  and 
u  =  36/80  as  shown  in  Figure  27.  The  slope  of  this  segment  is  given  by 


40(PC  -  Pb) 

b 


(43) 


Figure  25.  Resistance  reduction  coefficient  -  A  for  Method  SLAT  1 


d.  The  exponential  section  of  the  curve,  from  a  to  b  ,  is  of  the  form 


p  =  CuVn 


(44) 
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Figure  26.  Resistance  reduction  coefficient  -  B  for  Method  SLAT  1 

where  the  parameters  C,  n,  and  the  terminus  of  the  initial  linear  portion  pa  and 
ua  are  obtained  by  forcing  the  exponential  function  in  Equation  44  to  pass 
through  pb  and  ub  with  the  same  slope  s  as  segment  be  and  to  have  the  slope  kp  at 
the  terminus  of  the  initial  straight  line  segment  at  a.  This  results  in 


Pb 

n  =  - 

(45) 

3  Ub 

c-  Pt 

V" 

(46) 

1 

II 

a3 

(47) 

II 

e 

ft 

(48) 

(Note:  In  some  situations  Equations  45  through  48  may  result  in  unrealistic 
values  for  ua  and/or  pa.  If  this  occurs,  the  exponential  portion  is  omitted  and 
the  initial  linear  segment  is  extended  to  its  intersection  with  the  straight  line 
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Figure  27.  p-u  curves  by  Method  SLAT2 

section  be  or  until  the  maximum  resistance  pc  is  reached  whichever  comes  first. 
If  segments  Oa  and  be  do  not  intersect  at  realistic  values  of  pa  and  ua ,  segment 
be  is  omitted.) 


Murchison  and  O’Neill  (1984) 

Murchison  and  O'Neill  (1984)  simplified  the  process  of  Method  SLAT1  by 
replacing  the  three-part  p-u  curve  with  a  single  analytical  expression  as  follows. 


P_ 

Pu 


nA  tan  h 


\ 


kz 

”4 Pu ) 


(49) 


where 

pu  =  ultimate  lateral  soil  resistance  from  either  Equation  36  for  z  <  zcr  or 
Equation  37  for  z  >  zcr 

n  =  geometry  factor  =  1.5  for  tapered  piles  or  1.0  for  prismatic  piles 
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A  =  3-0.8 (z/b)  >  0.9  for  static  loads  or  =  0.9  for  cyclic  loading 
k  =  soil  stiffness  from  Table  2 
z  =  depth  at  which  the  p-u  curve  applies 
.  Several  illustrative  curves  for  this  method  are  shown  in  Figure  27. 


Synthesis  of  p-u  Curves  for  Piles  in  Clay 

Matlock  (1970) 

A  series  of  lateral  load  tests  on  instrumented  piles  in  clay  (Matlock  1970) 
were  used  to  produce  the  p-u  relationship  for  piles  in  soft  to  medium  clays 
subjected  to  static  lateral  loads  in  the  form 


/  N 1/3 

u 


\  UcJ 


(50) 


with  pu,  the  ultimate  lateral  resistance,  given  by  the  smaller  of 


Pu 


3 


b 


(51) 


for  a  wedge  failure  near  the  ground  surface,  or 

Pu  =  (52) 

for  flow  failure  at  depth;  and  uc ,  the  lateral  displacement  at  one-half  of  the 
ultimate  resistance,  given  by 

uc  =  2.5  e50b  (53) 

where 

y'  =  effective  unit  weight  of  the  soil 

su  -  shear  strength  of  the  soil 

J=  0.5  for  a  soft  clay  or  0.25  for  a  medium  clay 

e50  =  strain  at  50  percent  of  the  ultimate  strength  from  a  laboratory  stress- 
strain  curve 
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Typical  values  of  e50  are  given  in  Table  3.  The  depth  at  which  failure  transitions 
from  wedge  (Equation  51)  to  flow  (Equation  52)  is 


= 

cr  (y'b  +  J  su) 


(54) 


Figure  28.  p-u  curves  by  Method  CLAT1 

For  cyclic  loads,  the  basic  p-u  curve  for  static  loads  is  altered  as  shown  in 
Figure  28b.  The  exponential  curve  of  Equation  50  is  terminated  at  a  relative 
displacement  u/uc  =  3.0  at  which  the  resistance  diminishes  with  increasing 
displacement  for  z<zcr  or  remains  constant  for  z>zcr . 
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Table  3 

Representative  Values  of  e,0 

Shear  Strength  (psf) 

Percent 

250-500 

0.02 

500-1 ,000 

0.01 

1,000-2,000 

0.007 

2,000-4,000 

0.005 

4,000-8,000 

0.004 

Reese,  Cox,  and  Koop  (1975) 

Reese,  Cox,  and  Koop  (1975)  performed  lateral  load  tests  of  instrumented 
piles  to  develop  p-u  curves  for  piles  in  stiff  clay  below  the  water  table.  The  p-u 
curve  for  static  loading,  Figure  29,  consists  of  five  segments  determined  as 
follows. 


Figure  29.  p-u  curves  by  Method  CLAT2  for  static  loads 


a.  The  initial  linear  p-u  relationship  Oa  has  a  slope  equal  to  the  product  of 
soil  stiffness  k  (see  Table  4)  and  the  depth  z  to  the  location  at  which  the 
curve  applies. 
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Table  4 

Representative  Values  of  Lateral  Soil  Stiffness  k for  Piles  in  Clay 
for  Method  CLAT2 

Average  Undrained  Shear  Strength  (tsf)1  | 

Loading  Type 

0.5-1 

1  -2 

2-4 

Static  loading  -  ks  (pci) 

500 

1,000 

2,000 

Cyclic  loading  -  kr  (pci) 

200 

400 

800 

|| 1  Average  shear  strength  should  be  computed  from  the  unconsolidated  undrained  shear  strength  II 

II  of  the  soil  to  a  depth  of  five  pile  diameters. 

l| 

b.  The  second  segment  of  the  curve  is  parabolic  of  the  form 


P_  =  I 
Pc  2 


/  M/2 

u_ 

\Uc, 


(55) 


with  pc  taken  as  the  smaller  of 


Pc  = 


2  +  —  z  +  2.83  - 


s  b 

U 


(56) 


for  wedge  failure  near  the  ground  surface,  or 
Pc  =  11  Sub 


and 


uc  =  esob 


(57) 

(58) 


where  e50  =  strain  at  50  percent  of  ultimate  strength  from  a  laboratory  stress- 
strain  curve;  and  the  parameter  ,4  for  defining  pertinent  displacements  in 
Figure  29  is  obtained  from  the  curve  shown  in  Figure  30. 

c.  Points  a  and  b,  Figure  29,  are  joined  by  a  parabolic  curve  of  the  form 


P  _  1 


t  .  M/2 


\Uc, 


(59) 


where  pc  is  the  smaller  of 


Pc  = 


2  +  X.  ,  +  2.83  -L 

Su  SuJ 


s  b 

u 


(60) 
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for  a  wedge  failure  near  the  ground  surface,  or 

pc  =  (61) 

for  a  flow  failure  at  depth;  and 

uc  =  e50b  (62) 

where 

y'  =  effective  unit  weight  of  the  soil 
su  =  undrained  shear  of  the  soil 

e50  =  strain  at  50  percent  of  ultimate  strength  from  a  laboratory  stress-strain 
curve  (see  Table  3) 

The  parameter^  in  Figure  29  is  obtained  from  Figure  30. 
d.  Segment  cd  of  the  p-u  curve  in  Figure  29  is  of  the  form 


p_  =  I 

(  \ 

u 

1/2 

-  0.055 

/  \ 
“  -1 

Pc  2 

V  c) 

A  u 

\  ^ suc  y 

e.  A  second  linear  segment  joins  points  c  and  d  at  the  slope  indicated  in 
Figure  29  and  the  lateral  resistance  remains  constant  for  lateral  displace¬ 
ments  greater  than  u  =  18  Asuc . 
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The  p-u  curve  for  cyclic  loading  provided  by  Reese,  Cox,  and  Koop  (1975)  is 
illustrated  in  Figure  31.  The  curve  is  constructed  as  follows 

a.  The  initial  linear  p-u  relationship  Oa  has  a  slope  equal  to  the  product  of 
soil  stiffness  k  (see  Table  4)  and  the  depth  z  to  the  location  at  which  the 
curve  applies. 

b.  The  second  segment,  joining  points  a  and  b  (Figure  3 1)  is  an  exponential 
relationship  of  the  form 


-P-  =  1  - 

AC  Pc 


abs 


0.45m„ 


-  1 


2.5 


(64) 


where 

Ac  =  pressure  reduction  coefficient  from  Figure  30 


Figure  31 .  p-u  curve  by  Method  CLAT2  for  cyclic  loads 
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and 


«p  =  4.1  AcUc  ^ 

where  uc  is  given  by  Equation  62. 

c.  A  second  linear  p-u  relationship  joins  points  b  and  c  with  the  slope 

shown  in  Figure  31.  For  displacements  greater  than  u  =  \.%up,  the  lateral 
resistance  remains  constant. 


Reese  and  Welch  (1975) 


Reese  and  Welch  (1975)  performed  a  lateral  load  test  on  an  instrumented 
drilled  shaft  in  stiff  clay  above  the  water  table.  The  p-u  curve  obtained  from  the 
experimental  results  for  static  loads  is  shown  in  Figure  32.  The  curve  consists  of 
an  exponential  relationship  between  lateral  resistance  and  displacement  to  an 
ultimate  resistance,  after  which  the  resistance  remains  constant  for  further 
displacement.  The  requisite  exponential  relationship  is 


P_ 

Pu 


_1_ 

2 


(  \  1/4 

U 


\  Ucj 


(66) 


Figure  32.  p-u  curve  by  Method  CLAT3  for  static  loads 
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where 


pu  =  ultimate  resistance  obtained  as  the  smaller  from  Equation  51  with  J  -  0.5 
or  from  Equation  52 

uc  =  critical  lateral  displacement  obtained  from  Equation  88 

Th ep-u  curve  for  cyclic  loading,  shown  in  Figure  33,  is  constructed  as 
follows: 


a.  Values  of  p/pu  for  various  values  of  static  displacement  u/uc  are  com¬ 
puted  from  Equation  66. 

b.  The  displacement  for  cyclic  loading  for  each  value  of  p/pu  is  obtained 
from 


—  =  -i  +  9.6 

Uc  Uc 


\ 

P_ 

Puj 


log  (V) 


(67) 


LATERAL  DISPLACEMENT  -  u/uc 


Figure  33.  p-u  curve  by  Method  CLAT3  for  cyclic  loads 
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where 


ws  =  static  displacement  corresponding  to  p/pu 
N  =  number  of  cycles  of  load  application 


Reese  and  Sullivan  (1980) 

Each  method  for  p-u  curves  for  piles  in  clay  described  above  was  developed 
for  a  single  soil  profile;  hence  there  were  no  recommendations  provided  for 
transitioning  from  “soft”  clay  criteria  to  “stiff’  clay  criteria.  Sullivan  (1977)  and 
Sullivan,  Reese,  and  Fenske  (1979)  reexamined  the  data  for  soft  clays  (Matlock 
1970)  and  stiff  clays  (Reese,  Cox,  and  Koop  1975)  and  developed  a  unified 
criterion  (Reese  and  Sullivan  1980),  which  yields  computed  behaviors  that  are  in 
reasonable  agreement  with  both  soft  and  stiff  conditions.  However,  some 
judgement  on  the  part  of  the  user  is  required  in  selecting  appropriate  parameters 
for  use  in  the  prediction  equations. 

The  p-u  curve  by  the  unified  criteria  for  static  loading,  illustrated  in 
Figure  34,  consists  of  an  initial  linear  segment  0a,  an  exponential  segment  ab,  a 
second  linear  segment  be  and  a  constant  lateral  resistance  for  large  displace¬ 
ments.  The  curve  for  static  loading  at  a  particular  depth  z  is  constructed  as 
follows: 
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(68) 

(69) 


where 

av  =  average  effective  vertical  stress  over  the  depth  z 
ca  =  average  cohesion  over  the  depth  z 
c  =  cohesion  at  depth  z 
b  =  pile  diameter 

(2)  For  z  >  \2b ,  the  ultimate  resistance  is 

pu  =  12 cb  (70) 

b.  Compare  the  properties  of  the  soil  profile  under  analysis  with  those  listed 
in  Table  5  and  select  the  values  of  parameters  Aand  F  to  be  used  in  the 
following  calculations. 

c.  The  p-u  relationship  for  the  initial  linear  segment  is 


—  =  kz 
Pu 

where  it  is  a  stiffness  parameter  from  Table  6  (see  also  Table  4). 

d.  The  exponential  segment  ab  is  obtained  from 

/  \  1/2 

JL  =  0.5  JL 

Pu  l  Uc) 
with 

uc  =  Aesob 


(71) 


(72) 


(73) 
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Table  5 

Curve  Parameters  for  Method  CLAT4 


Clay  Description 


Table  6 

Representative  Values  of  Jr  for  Method  CLAT4 


Cohesion  (os 


200-500 


500-1,000 


1 ,000-2,000 


2,000-4,000 


4,000-8,000 


p 


(74) 


Pu 


=  F  +  (1 


for  z  <  12 b  .  For  z  >  12 b  ,  p/pu  -  1. 


The  p-u  curve  by  the  unified  method  for  cyclic  loading,  Figure  35,  also 
consists  of  an  initial  linear  segment,  followed  by  an  exponential,  variation  of  p 
with  u,  a  second  linear  segment,  and  a  constant  resistance  for  large  displace¬ 
ments.  Construction  of  the  curve  for  cyclic  loading  follows  the  same  steps  as  for 
the  static  curve,  with  the  exceptions  that  the  exponential  segment  terminates  at  a 
resistance  equal  to  one  half  of  pu ,  the  second  linear  segment  terminates  at  a 
displacement  u  =  20uc ,  and  the  constant  resistance  for  u  >  20 uQ  is  given  by 


(75) 


for  z  <  12 b .  For  z  >  12 b  ,  p/pu  =  1. 


Figure  35.  p-u  curve  by  Method  CLAT4  for  cyclic  loading 


O’Neill  and  Gazioglu  (1984) 

Although  the  procedure  presented  as  Method  CLAT4  attempted  to  provide  a 
unified  criterion  for  all  clays,  the  procedure  requires  the  user  to  select  parame¬ 
ters  a  priori  which  essentially  convert  the  method  to  a  soft-clay-like  method  or  to 
a  stiff-clay-like  method  (O'Neill  and  Gazioglu  1984).  O’Neill  and  Gazioglu 
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reexamined  the  data  utilized  in  developing  the  previous  methods  as  well  as  the 
results  of  other  tests  of  instrumented  laterally  loaded  piles  in  clay  to  produce  an 
integrated  procedure  for  p-u  curves  for  piles  in  clay.  The  method  attempts  to 
incorporate  continuum  effects  and  relative  pile/soil  stiffness  characteristics 
which  were  not  explicitly  accounted  for  in  the  previous  procedures. 

O’Neill  and  Gazioglu,  as  well  as  other  researchers,  reasoned  that  there  exists 
a  critical  length  of  pile  such  that  longer  piles  no  longer  influence  the  pile  head 
behavior.  This  critical  length  is  presented  as 


El 


Esb 


1/2 


0.286 


(76) 


where  EI=  flexural  stiffness  of  the  pile;  and  Es  =  “perhaps  a  secant  Young's 
modulus  at  a  deviator  stress  level  of  one-half  of  the  deviator  stress  at  failure  in 
undrained  triaxial  compression,”  with  all  length  units  in  inches.  Es  is  evaluated 
for  the  average  UU  triaxial  shear  strength  between  the  ground  surface  and  the 
depth  Lc .  Hence  the  determination  of  Lc  is  an  iterative  process.  Representative 
values  of  Es  are  given  in  Table  7. 


Table  7 

Soil  Modulus  for  Method  CLAT5 

Undrained  Shear  Strength  c(psf) 

Soil  Modulus  (psi) 

<500 

50 

500  -  1 ,000 

50-150 

1,000-2,000 

150-450 

2,000  -  4,000 

450-1,500 

4,000  -  8,000 

1,500-5,000 

>8,000 

5,000 

The  reference  lateral  displacement  uc  (see  Equations  53,  58,  62,  and  73)  is 
obtained  from 


uc  =  A'e50b 


1/2 


'e? 


1/8 


(77) 


where  A'  =  constant  taken  as  0.8  by  O’Neill  and  Gazioglu;  Es  =  soil  modulus 
from  Table  7  for  the  depth  of  interest;  and  all  length  units  are  in  inches. 

The  ultimate  soil  resistance  is  expressed  as 

pu  =  FNpcb  (78) 
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where  F  =  reduction  factor  from  Table  8  for  the  appropriate  loading  condition; 
and  Np  is  given  by 


N  =  3+  6—  <9  (79) 

p 


Table  8 

Soil  Degradability  Factors 

UU  Triaxial  Compression  Failure  Strain  j| 

Factor 

Loading  Condition 

<0.02 

0.02-0.06 

>0.06 

F 

Static 

0.50 

0.75 

1.00 

F 

Cyclic 

0.33 

0.67 

1.00 

where  zcr  -  LJ4,  indicating  a  transition  from  wedge  type  failure  to  flow  failure  at 
a  depth  equal  to  one  fourth  the  critical  length  which,  unlike  the  previous  meth¬ 
ods,  reflects  the  relative  pile/soil  stiffness. 

The p-u  curves  for  static  load  (Figure  36)  and  cyclic  loads  (Figure  37)  have 
an  initial  exponential  relationship  between  lateral  resistance  and  displacement 
given  by 


0.387 


V 


u 


c) 


(80) 


with  the  static  curve  terminating  at  p/pu  =  1  and  the  cyclic  curve  at  p/pu  -  0.5. 
The  second  linear  section  of  the  curves  terminates  at 


-F,)f 

Pu  Zcr 

and  the  cyclic  curve  at 


where  Fs  and  Fc  are  given  in  Table  8. 


(81) 


(82) 


Algorithm  for  Analysis  of  Laterally  Loaded  Piles 

The  p-u  curves  described  in  the  preceding  paragraphs  were  derived  on  the 
assumption  that  the  lateral  resistance  p  at  any  point  on  the  pile  is  a  function  only 
of  the  lateral  displacement  u  at  that  point  (i.e.,  the  Winkler  assumption).  For  this 
assumption  and  the  one-dimensional  model  of  the  pile-soil  system  shown  in 
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Lateral  Displacement  -  u/uc 


Figure  36.  p-u  curve  by  Method  CLAT5  for  static  loading 


Figure  37.  p-u  curve  by  Method  CLAT5  for  cyclic  loading 


Figure  38,  the  governing  differential  equation  for  bending  in  the  x-z  plane  of  a 
prismatic,  linearly  elastic  pile  is 
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d2 

'  El  fa' 

+  d_ 

dz 2 

dz*  ) 

dz 

dz 

-  p(z,u)  =  0 


(83) 


where 

E  =  modulus  of  elasticity  of  the  pile  material 

/  =  moment  of  inertia  of  pile  cross  section  about  an  axis  perpendicular  to 
the  x-z  plane 

P(z)  =  axial  compressive  force  in  the  pile  at  z 

p(z,u)  =  lateral  resistance  which  is  a  function  of  both  position  z  on  the  pile  and 
the  lateral  displacement  u  at  z 
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Because  the  displacement  u  must  be  known  before  the  lateral  resistance  p  can 
be  evaluated,  numerical  iterative  solutions  of  Equation  83  are  required.  The 
most  common  approach  is  to  represent  the  pile-soil  system  by  a  discretized 
model  such  as  illustrated  in  Figure  38  where  the  displacements  and  forces  are 
evaluated  at  a  finite  number  of  points  (Matlock  and  Reese  1962,  Dawkins  1982). 
The  solution  proceeds  as  a  succession  of  trials  and  corrections  until  forces  and 
displacements  are  compatible  at  every  node. 


Observations  of  System  Behavior 

Obtaining  the  numerical  solutions  described  above  is  expedited  by  replacing 
the  nonlinear  p-u  curve  by  equivalent  linearly  elastic  springs  during  each 
iteration.  The  stiffnesses  of  these  linear  springs  are  evaluated  as  the  secant  to 
p-u  curve  for  the  displacements  calculated  during  the  preceding  iteration.  It  is  to 
be  noted  that  the  ultimate  lateral  resistance  tends  to  increase  with  depth,  while 
the  pile  displacements  decrease  as  z  increases.  Hence  it  can  be  concluded  that 
the  secant  stiffness  of  the  lateral  resistance  increases  with  depth  below  the  pile 
head.  Consequently  if  the  distribution  of  soil  stiffnesses  along  the  pile  can  be 
determined,  the  behavior  may  be  evaluated  for  working  loads  without  the  need 
for  iterative  solutions. 

The  second  term  in  Equation  83  represents  the  interaction  of  the  axial  load  in 
the  pile  with  the  lateral  displacement  to  increase  the  bending  moments  in  the  pile 
(the  “beam-column”  effect).  Distribution  of  the  axial  load  is  not  influenced 
significantly  by  lateral  loading;  hence,  the  distribution  may  be  determined  using 
the  axial  load  analysis  techniques  described  in  Chapter  2  before  a  lateral  load 
analysis  is  performed.  In  usual  pile-soil  systems,  the  beam-column  effect  is 
small  and  conservative  estimates  of  its  influence  may  be  obtained  by  taking  the 
axial  force  in  the  pile  equal  to  the  applied  head  load  P0. 


Linearly  Elastic  Analyses 

A  linear  relationship  between  lateral  resistance  and  lateral  displacement  is 
expressed  by 

p(z,u)  =  -Ez(z)u  (84) 

where  E(z)  =  secant  stiffness  of  the  lateral  resistance  and  the  minus  sign  indi¬ 
cates  that  the  resistance  is  opposite  in  direction  to  the  displacement.  For  a 
prismatic  pile  with  constant  axial  force  and  a  linear  resistance  function,  Equation 
83  becomes 


£u 
dz 2 


+  Ez(z)u  =  0 


(85) 
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Even  with  the  simplifications  inherent  in  Equation  85,  explicit  closed  form 
solutions  are  not  possible  unless  the  lateral  resistance  stiffness  is  constant  with 
depth.  Although  the  lateral  resistance  is  not  constant,  insight  into  the  lateral  load 
behavior  can  be  gained  from  the  solution  for  constant  Ez.  In  this  case  the 
differential  equation  becomes 


El  —  +  P 


+  Eu 


The  solution  of  Equation  86  (Hetenyi  1946)  is 

u  =  e  -az,T  [  c ,  cos—  +  C2sin— )  +  e^T  {  C3 cos—  +  C4sin— ]  (87) 

1  T  T)  T  T) 

where  C, ,  C2 ,  C3 ,  and  C  are  constants  to  be  determined  from  conditions  at  the 
pile  head  and  tip;  and 


r4  _  4 m 

(88) 
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II 
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<3 

(89) 

b2  =  1  +  X2 

(90) 

( p  r2\ 

(91) 
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For  a  sufficiently  long  pile  (see  Page  51)  conditions  at  the  pile  tip  have  no 
effect  on  the  response  at  the  pile  head.  In  this  case  the  relationship  between 
loads  (M0 ,  V0)  and  displacements  (u0 ,  p0)  at  the  pile  head  may  be  expressed  as 


For  a  value  of  A.2  =  1/2  ,  the  determinant  of  the  coefficient  matrix  in  Equation  92 
is  zero,  indicating  that  an  axial  head  load 

Po  =  ( Efil)m  (93> 

is  the  buckling  load  for  the  long  pile.  Hetenyi  (1946)  gives  the  buckling  load  for 
a  finite  length  pile  with  both  ends  fixed  against  displacement  as  less  than  or 
equal  to 
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(94) 


f.  2(£IH)'n 

L 

For  usual  pile-soil  systems  the  axial  head  load  must  be  significantly  less  than  the 
value  given  by  Equation  93  in  order  to  prevent  overstress  of  the  pile  material  at 
the  head.  Hence  X2  is  always  considerably  less  than  1,  so  that  the  beam-column 
effect  may  be  neglected  and  the  second  term  of  Equations  83,  85,  and  86  may  be 
omitted.  Consequently,  the  governing  differential  equation  for  elastic  analyses 
becomes 

El  —  +  E(z)u  =  0  (95) 

dz 4 


Variation  of  Lateral  Resistance  Stiffness 

In  order  to  utilize  Equation  95  for  analysis  of  laterally  loaded  piles  the 
variation  and  magnitude  of  the  lateral  resistance  stiffness  must  be  known. 
Terzaghi  (1955)  provides  “coefficients  of  horizontal  subgrade  reaction”  for 
constant  E.  for  clay  (Table  9)  and  for  Es  varying  linearly  with  depth  (E/z)  =  n^z) 
for  sand  (Table  10).  The  work  of  Skempton  (1951)  has  been  extended  to  piles  in 
homogeneous  soft  clays  to  evaluate  Es  as 


E.  = 


5es0b 


(96) 


Table  9 

Values  of  £,for  1-ft-Wide  Piles  in  Precompressed  Clay 

Consistency  of  Clay 

Stiff 

Very  Stiff 

Hard 

Values  of  cu  (tsf) 

1-2 

2-4 

>4 

Range  of  E$  (tcf) 

33-67 

67-133 

>133 

Proposed  value  of  E ,  (tcf) 

50 

100 

200 

Table  10 

Values  of  Constant  of  Horizontal  Subgrade  Reaction  nh  for  a  1-ft- 
Wide  Pile  in  Sand 

Relative  Density 

Loose 

Medium 

Dense 

nh  (tcf)  for  dry  or  moist  sand 

7 

21 

56 

nh  (tcf)  for  submerged  sand 

4 

14 

34 
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The  p-u  curves  discussed  previously  provide  an  indication  of  the  variation  of 
Es  with  depth.  For  both  clays  and  sands,  the  ultimate  resistance  used  for  con¬ 
structing  the  p-u  curves  increases  with  depth  near  the  ground  surface.  At  some 
depth,  the  ultimate  resistance  for  clays  reaches  a  limiting  value.  Considering 
that  the  lateral  displacements  of  the  pile  decrease  exponentially  with  depth  (see 
Equation  87),  it  seems  probable  that  the  secant  stiffness  of  the  lateral  reaction  for 
both  sands  and  clays  increases  as  some  exponential  function  with  depth  of  the 
form  (Matlock  and  Reese  1962). 

Ez(z)  =  K0  +  Kz"  (97) 

The  stiffness  of  sand  at  the  ground  surface  will  be  zero;  hence,  K0  will  be  zero. 
Although  the  p-u  curves  for  clays  suggest  that  K0  will  not  be  zero,  it  is  conserva¬ 
tive  to  take  K0  as  zero  for  these  materials.  Therefore,  the  elastic  solutions 
presented  in  the  following  paragraphs  will  take  K0  to  be  zero  in  all  cases. 

Although  explicit  solutions  of  Equation  95  exist  for  Es  constant  with  depth 
(k  =  0  ,  Es  =  K  in  Equation  97),  all  solutions  are  presented  in  graphical  form. 
Following  the  procedures  of  Matlock  and  Reese  (1962)  the  following  non- 
dimensional  parameters  are  defined 


m 

K 


(98) 


Z  = 


z_ 

T 


L 

T 


(99) 

(100) 


u  = 


VJ_ 

El 


AU{Z) 


MT 2 

irw 


(101) 


where  V0  and  M0  are  the  applied  head  shear  and  moment,  respectively;  and  Au( Z) 
and  BU(Z)  are  nondimensional  functions  of  the  nondimensional  depth  Z. 

Substitution  of  Equations  97  through  101  into  Equation  95  yields 
dU  (Z) 

- u—  +  ZnA(Z)  =  0  (102) 

dZ 4 

and 

d4B  (Z) 

- ^  +  ZnBu(Z)  =  0  (103) 

dt 
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With  solutions  of  Equations  102  and  103  available,  the  displacement  at  any 
point  on  the  pile  is  obtained  from  Equation  136  and  the  slopes  of  the  pile,  the 
bending  moments,  and  shears  are  evaluated  from  the  derivatives  of  u  as 


’2  MT 


,  v  T1 
p  =  —  =  AS(Z)  +  ^  BS(Z) 
dz  El  s  .  El  1 


(104) 


rjAJZ)  * mjijz) 

dr 


(105) 


d*u  A/ 

V  =  EI  =  V/lv(Z)  +  BV(Z) 

dz  T 


(106) 


The  various  functions  of  A  and  B  are  plotted  in  Appendix  A  (Figures  B1  through 
B24)  for  constant,  linear,  and  parabolic  variations  of  E,  with  depth  («  =  0,  1,  and 
2,  respectively,  in  Equation  96). 


Pile  Head  Stiffness  Coefficients 
for  Lateral  Loading 

Pile  head  stiffness  coefficients  for  lateral  loading  are  obtained  by  inverting 
Equations  100  and  101  with  Z=  0  as  shown  in  the  following  matrix  equation 


EI  ..... 

'  K 

~TBuo 

4 

Uo 

Mo. 

r3M„  - 

o 

S 

l 

T2Auo 

P°. 

(107) 


where 


AU0=AU(Z  =  0) 


BU0  =  BU(Z  =  0) 


Aso=A/Z  =  0) 


Bso  =  B/Z  =  0) 


The  coefficients  Au0 ,  Buo ,  Aso ,  and  Bso  are  shown  for  various  relative  pile  lengths 
Zmax  in  Appendix  B  (Figures  B25  through  B28).  The  following  items  should  be 
noted:  (a)  As0  =  Buo  ;  (b)  piles  with  Zmax  <  2  may  be  treated  as  rigid  (see  page 
76);  and  (c)  The  A  and  B  coefficients  remain  constant  for  Zmax  >  4. 
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As  noted  on  page  59,  piles  with  Zmax  >  4  may  be  treated  as  inflexible.  In  this 
case  the  lateral  displacement  at  any  depth  may  be  expressed  in  terms  of  the  pile 
head  displacements  as 

u  =  u0  +  pjz  (108) 

and  the  soil  resistance  at  that  point  is 

p  =  Kznu  =  Kzn{u0  +  poz)  (109) 

From  an  equilibrium  analysis  of  the  rigid  pile,  the  head  shear  and  moment  are 
given  in  terms  of  the  head  displacements  by 


K 

•  =  KLnA 
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Evaluation  of  Linear  Lateral  Soil  Resistance 

In  order  to  apply  the  linearized  solutions  described  in  the  preceding  para¬ 
graphs,  the  variation  and  magnitude  of  the  lateral  soil  resistance  stiffness  must 
be  evaluated.  Terzaghi  (1955)  provides  estimates  of  clay  soil  stiffness  constant 
with  depth  ( n  =  0)  and  sand  soil  stiffness  varying  linearly  with  depth  (n  =  1)  as 
shown  in  Tables  9  and  10.  (Note:  Terzaghi  states  that  the  soil  stiffness  values 
are  for  a  “1  foot  wide  pile”  and  in  order  to  apply  these  values  to  piles  of  different 
widths  the  stiffness  for  the  1-ft-wide  pile  must  be  divided  by  the  actual  width  of 
the  pile.  In  order  to  utilize  the  resulting  “horizontal  subgrade  modulus”  in  the 
linearized  analysis,  Terzaghi's  modulus  must  be  multiplied  by  the  width  of  the 
pile  in  contact  with  the  soil  (see  Hetenyi  1941).  Consequently  the  moduli  given 
by  Terzaghi  may  be  used  without  alteration  as  the  value  of  K  in  the  linearized 
equations.) 

Because  the  laterally  loaded  pile-soil  system  is  highly  nonlinear,  particularly 
under  large  loads,  immutable  pile  head  stiffness  coefficients  do  not  exist. 
Although  the  soil  stiffness  moduli  given  in  Tables  9  and  10  can  be  used  to 
evaluate  explicit  coefficients,  these  values  must  be  interpreted  as  only  first 
approximations.  Higher  approximations  may  be  obtained  by  combining  the 
nondimensional  solutions  with  the  nonlinear  p-u  curves  discussed  earlier  (Reese, 
Cooley,  and  Radhakrishnan  1984,  “Executive  Summary....”)  as  outlined  in  the 
following  steps. 
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a.  Evaluate  p-u  curves  for  the  appropriate  soil  profile.  These  curves  should 
be  closely  spaced  in  the  top  10  to  20  pile  diameters. 

b.  Estimate  a  variation  and  lateral  stiffness  (i.e.,  K  and  n  in  Equation  97)  for 
the  soil  profile  using  Terzaghi's  soil  moduli.  (Reese,  Cooley,  and 
Radhakrishnan  (1983)  suggest  that  a  value  of  T  (Equation  98)  be 
assumed;  Terzaghi's  moduli  provide  a  means  for  this  assumption.) 

c.  Evaluate  the  deflections  at  the  locations  of  the  p-u  curves  in  step  a  using 
the  appropriate  nondimensional  curves  for  head  loads  in  the  working 
range. 

d.  Determine  the  slope  of  a  secant  line  from  the  p-u  curve  for  the  deflection 
calculated  for  each  location.  This  establishes  the  soil  modulus  E.  at  each 
p-u  curve  location  and  allows  E,  to  be  plotted  versus  depth  z. 

e.  Revise  the  variation  and  lateral  soil  stiffness  (i.e.,  new  K  and  n  in  Equa¬ 
tion  97)  to  best  approximate  the  curve  of  E.  versus  z  in  step  d. 

f.  Repeat  steps  c,  d,  and  e  until  convergence  is  achieved. 

g.  Use  the  final  values  of  K  and  n  to  calculate  the  pile  head  stiffness  coeffi¬ 
cients  in  Equations  107  and  1 10. 
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4  Algorithm  for  Analysis  of 
Torsionally  Loaded  Single 
Piles 


Three-dimensional  analysis  of  a  single  pile  requires  a  relationship  between 
the  resistance  of  the  soil  and  the  torsional  displacement  of  the  pile.  There  has 
been  only  limited  investigation  (O'Neill  1964,  Poulos  1975,  Scott  1981,  Stoll 
1972)  of  this  torque-twist  relation  because  its  effect  is  small  compared  to  the 
axial  and  lateral  effects.  Until  more  detailed  data  are  available,  the  following 
simplistic  relationship  should  be  used. 

It  is  assumed  that  the  soil  is  a  radially  linearly  elastic,  homogenous  medium, 
that  the  pile  is  prismatic  and  linearly  elastic,  and  that  the  resistance  of  the  soil  at 
any  point  is  a  function  only  of  the  torsional  displacement  of  the  pile-soil  inter¬ 
face  at  that  point.  Under  these  assumptions  the  soil  in  any  plane  perpendicular 
to  the  axis  of  the  pile  is  in  a  state  of  plane,  pure  shear.  The  theory  of  elasticity 
solution  for  this  case  yields 


T 


(111) 


and 


V?2 

2  G/3 


(112) 


where 

t  =  shear  stress  at  a  radial  distance  r  from  the  centerline  of  the  pile 
r0  =  shear  stress  at  the  pile-soil  interface 
R  =  radius  of  the  pile 
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V  =  displacement  perpendicular  to  the  radial  direction  at  r 
Gs  =  shear  modulus  of  the  soil 


If  there  is  no  slippage  between  the  pile  and  soil  at  the  interface,  the  tangential 
displacement  of  a  point  on  the  interface  is 

v  =  i?0  (113) 

where  6  is  the  rotation  of  the  pile.  And,  finally,  the  required  relation  is 


T 

-2-  =  2G 

e 


(114) 


The  linear  relationship  between  surface  shear  and  pile  rotation  represented  by 
Equation  1 13  is  assumed  to  terminate  when  the  surface  shear  r0  reaches  a  limit 


of 


t  =  k  o  tan  6 

ou  o  v 


for  sands,  or 


(115) 


(116) 


for  clays,  where 

vou  -  ultimate  surface  shear  resisting  rotation  of  the  pile  about  its  longitudinal 
axis 

k0  =  at-rest  pressure  coefficient 
6  =  angle  of  pile-soil  interface  friction  for  sand 
o0  =  vertical  effective  stress 

a-  an  adhesion  factor  which  may  be  obtained  from  Figure  12 

su  =  shear  strength  of  clay.  The  resistance  to  rotation  remains  constant  at  vou 
for  additional  rotational  displacement  as  shown  in  Figure  39 


Elastic  Analysis 

So  long  as  the  surface  shear  is  less  than  xou ,  the  entire  pile-soil  system  is 
linearly  elastic.  The  governing  differential  equation  for  torsional  response  of  the 
linear  pile-soil  system  is 
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Figure  39.  Proposed  torsional  shear  -  rotation  curve 


GJ  -  4 tiR2G  (z)Q  =  o 

dz2 

where 


(117) 


G  =  shear  modulus  of  the  pile  material 

j  —  torsional  area  property  of  the  pile  cross  section  (polar  moment  of  inertia 
for  a  circular  section) 

Because  Equation  1 16  is  identical  in  form  to  the  differential  equation  for  an 
axially  loaded  pile,  pile  head  torque-twist  stiffness  may  be  obtained  from  the 
equations  and  procedures  appearing  on  pages  A3-A6  of  Appendix  A  by  perform¬ 
ing  the  following  substitutions: 

a.  In  Equations  A4  and  A5  (in  Appendix  A),  define  T,  =  (GJ/4tcR2Gs)112  for 
Gs  constant  with  depth. 

b.  In  Equations  A7  and  A1 1,  replace  EA  with  GJ ;  replace  w0  with  60  (the 
twist  angle  at  the  pile  head);  and,  replace  P0  with  M0  (the  torsional 
moment  at  the  pile  head). 
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c. 


In  Equation  A16,  replace  EA  with  GJ ;  replace  w(z)  with  0(z);  and, 
replace  P0  with  M0. 


d.  In  Equation  A18,  replace  EA  with  GJ ;  and,  define  Kf  such  that  Gs(z )  - 
Kfzn  for  Gs  varying  with  depth. 

(Note:  Scott  (1981)  indicates  that  the  torsional  resistance  to  twist  at  the  pile  tip 
may  be  included  as  was  done  for  tip  reaction  for  the  axially  loaded  pile.  How¬ 
ever,  in  most  situations  the  tip  resistance  against  twist  will  be  negligible.) 
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Three-Dimensional  System 

Figure  40  illustrates  the  coordinate  system,  forces,  and  displacements  at  the 
pile  head  which  must  be  considered  in  a  three-dimensional  analysis.  The  x-  and 
y-axes  are  the  principal  axes  of  the  pile  cross  section  and  the  z-axis  is  the 
longitudinal  axis  of  the  pile.  Forces  and  displacements  are  assumed  to  have 
positive  senses  in  the  positive  coordinate  directions  (“right-hand  rule  for 
moments  and  rotations). 

For  a  linearly  elastic  system,  the  forces  and  displacements  are  related  by 


The  b  coefficient  matrix  array  is  the  pile  head  stiffness  matrix  and  the  individual 
elements  by  are  obtained  from  Equations  A21, 107,  and  108. 
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where 


E  =  modulus  of  elasticity  of  pile  material 

I  =  moment  of  inertia  of  pile  cross  section  about  y-axis 

Ta  =  length  parameter  for  lateral  loading  in  the  x-z  plane 

Equation  98  coefficients  Auo ,  As0 ,  Bu0 ,  and  Bso  are  obtained  from  Figures  B25 
through  B28  with  Zmax  =  L/Tn ;  terms  appearing  in  Equation  122  are  defined  in 
Chapter  2;  and,  terms  in  Equation  123  are  defined  in  Chapter  4.  The  remaining 
elements  of  the  pile  head  stiffness  matrix,  b22 ,  b2A ,  b55 ,  and  bu  ,  are  evaluated 
for  bending  in  the  y-z  plane. 


Pile  Head  Fixity 

If  the  pile  head  is  attached  to  the  supported  structure  so  that  the  displace¬ 
ments  of  the  pile  head  and  the  point  of  attachment  on  the  structure  undergo 
identical  displacements,  the  stiffness  matrix  as  shown  in  Equation  118  may  be 
included  as  a  part  of  the  overall  system  stiffness  without  alteration.  In  most 
installations,  the  pile  head  and  the  supported  structure  will  experience  the  same 
translational  displacements  ( u,v,w ).  However,  the  method  of  connection  may 
permit  relative  rotation  between  the  structure  and  the  pile.  To  illustrate  the 
effect  of  relative  rotation  of  the  pile  and  structure,  the  two-dimensional  system 
shown  in  Figure  38  is  used.  The  relationship  between  the  head  forces  and 
displacements  is 


u 


O 


Po 


VT3 

-2—  A  + 
El 


VT 2 

-2—  A,„  + 
El  s° 


MT2 

0  B 

El 

(123) 

MT 

0  B 

El  s° 

(124) 

The  attendant  inverse  relationship,  considering  only  the  terms  associated  with 
lateral  loading  and  the  notation  of  Equation  1 19, 
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Pinned-Head  Pile 


If  the  pile-to-structure  connection  is  such  that  no  moment  is  transmitted 
through  the  connection,  then  M0  will  be  zero.  For  a  unit  lateral  translation, 
u0  =  1,  Equations  123,  124,  and  125  yield 

V  T3 

1  =  A*°  (126) 

or 

h,  =  a  j_ 

11  '  t3  A  (127) 

*  UO 

and,  the  resulting  rotation  of  the  pile  head, 


P0  =  ~ 


1  As 


T  A., 


(128) 


with  b[5-  b5l  =  b55  =  0  .  (The  prime  superscript  denotes  the  pinned  head 
condition.) 


Partial  Fixity  at  Pile  Head 


Frequently  the  pile-structure  connection  permits  a  limited  relative  rotation 
before  moment  resistance  at  the  pile  head  is  developed.  To  simulate  the  partial 
fixity,  it  is  assumed  that  moment  resistance  develops  at  a  reduced  rate  propor¬ 
tional  to  the  degree  of  fixity  f(o<f<  1) .  To  evaluate  the  stiffness  elements  bn 
and  b$l  for  partial  fixity,  a  unit  value  of  u0  is  imposed  at  the  pile  head  and  the 
rotation  is  allowed  to  increase  to  (1  -/)  p'0.  Similarly,  to  evaluate  bI5  and  b5S , 
u0  =  0  and  p0  =  fp'0  are  imposed.  The  resulting  stiffness  elements  are 
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Free-Standing  Pile  Segment 


A  portion  of  the  pile  may  extend  above  the  ground  surface  as  illustrated  in 
Figure  41  for  bending  in  the  x-z  plane.  Although  the  free-standing  segment 
may  be  considered  as  a  part  of  the  structure  supported  by  the  pile,  it  may  be 
advantageous  to  combine  the  free-standing  and  embedded  segments  to  express 
the  relationship  between  forces  and  displacements  at  the  pile  head,  point  o  in 
Figure  41,  and  eliminate  consideration  of  the  effects  at  the  ground  surface. 


Figure  41 .  Linearly  elastic  pile/soil  system  with  free-standing  segment 


The  free-standing  segment  is  not  subjected  to  the  effects  of  surrounding  soil; 
hence,  the  relationship  between  end  forces  and  end  displacements  for  the  free¬ 
standing  portion  is  required.  Using  conventional  beam-column  element  theory 
for  the  two-dimensional  system  in  Figure  41,  the  free-standing  segment’s 
STIFFNESS  matrix  can  be  written  as: 
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where  E  ,A  ,Iy,  and  l  are  the  modulus  of  elasticity,  cross-sectional  area, 
moment  of  inertia  and  length,  respectively,  of  the  free-standing  segment,  and  / 
is  the  degree  of  fixity  between  the  pile  head  and  structure  at  point  o  .  For 
convenience,  the  matrices  in  Equation  1 32  are  partitioned  as  indicated  by  the 
dashed  lines  and  this  equation  may  be  written  symbolically  as 


* 

b~  b~ 

u~ 

oo  "os 

0 

< 

b~  b~ 

u~ 

so  ss 

Is) 

(133) 


Similarly,  the  relationship  between  forces  and  displacements  at  point  5  (Fig¬ 
ure  41)  may  be  written  as  (see  Equation  118) 
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Because  there  are  no  external  forces  at  point  s.  Equations  133  and  135  may  be 
combined  to  obtain 


=  K0u~o 


(136) 


P.-o-  b-v. :  *  (b-  *  b;)u ; 

From  Equation  137 

u~s  =  -(Ks +  b;T  Kou~o 


(137) 


(138) 


And  substitution  of  Equation  138  into  Equation  136  yields 


F0  = 


u;  =  [b~]u~0 


(139) 


where 


[B~] 


bn  bg  0 


b$$  0 


0  0  633' 


(140) 


is  the  pile  head  stiffness  matrix  for  bending  in  the  xz-plane.  A  similar  operation 
is  required  for  bending  in  the  yz-plane.  The  torsion  stiffness  coefficient  is  given 
by 


b 


t  /  f 
66 


+ 


-1 


(141) 


where  b66  is  the  torsional  coefficient  for  the  embedded  segment  from  Chapter  4, 
and  G  and  J  are  shear  modulus  and  torsional  area  moment  of  inertia,  respec¬ 
tively,  for  the  free-standing  segment. 

Alternatives  for  Evaluating  Pile  Head  Stiffnesses 

The  most  reliable  means  of  evaluating  the  pile  stiffness  is  from  field  tests  of 
prototype  piles.  Although  the  coefficients  relating  lateral  head  loads  and 
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displacements  may  be  evaluated  from  lateral  load  tests,  such  tests  are  not 
routinely  performed. 

For  complex  soil  conditions  and/or  nonprismatic  piles  which  are  not  readily 
approximated  by  one  of  the  procedures  for  linearly  elastic  systems  discussed 
previously,  the  pile  head  stiffness  matrix  may  be  obtained  with  the  aid  of 
computer  programs  such  as  CBEAMC,  CAXPILE,  or.COM624. 
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6  Analysis  of  Pile  Groups 


Although  isolated  single  piles  may  be  encountered  in  some  applications,  it  is 
more  common  that  a  structure  foundation  will  consist  of  several  closely  spaced 
piles  (many  building  codes  require  a  minimum  of  three  piles  in  a  group).  The 
structure/pile/soil  system  is  highly  indeterminate  and  nonlinear.  Historically, 
design  methods  have  been  based  on  numerous  simplifying  assumptions  that 
render  the  analytical  effort  tractable  for  hand  computations.  The  advent  of  the 
computer  has  allowed  solutions  to  be  obtained  in  which  many  of  the  simplifica¬ 
tions  of  the  classical  design  methods  are  no  longer  necessary.  Synopses  of  some 
of  the  classical  methods  and  more  complete  descriptions  of  the  computer-based 
techniques  are  presented  below. 


Classical  Methods  for  Pile  Group  Analysis 

All  of  the  classical  methods  assume  that  the  pile  cap  (or  super-structure)  is 
rigid  and  that  all  loads  are  resisted  only  by  axial  forces  in  the  piles.  These 
methods  attempt  to  allocate  the  superstructure  loads  to  individual  piles  through 
the  equations  of  static  equilibrium.  No  direct  attempt  is  made  to  determine  the 
deformations  of  the  system. 


Moment-of-lnertia  (Simplified 
Elastic  Center)  Method 

A  complete  description  of  the  Elastic  Center  method  is  given  by  Andersen 
(1956).  For  the  simplified  procedure  presented  here,  it  is  assumed,  in  addition  to 
a  rigid  cap,  that  only  vertical  loads  are  applied  to  the  cap,  that  all  piles  are  verti¬ 
cal,  that  all  piles  have  the  same  axial  stiffness  (EA/L),  and  that  the  magnitudes 
of  the  axial  loads  in  the  piles  vary  linearly  with  distance  from  the  centroid  of  the 
pile  group.  The  axial  load  at  the  head  of  the  i*  pile  is  given  by 


MyXj  ±  Mxyt 


(142) 
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where 


V  =  resultant  vertical  load  on  the  cap 
n  =  number  of  piles  in  the  group 

Ix,Iy  =  moments  of  inertia  about  x-  andy-axes,  respectively,  through  the 
centroid  of  the  piles  which  are  treated  as  point  (unit)  areas 

Mx,  My  =  moments  of  the  vertical  loads  on  the  cap  about  the  x-  and  y-axes, 
respectively 


Culmann's  Method 

The  method  attributed  to  Culmann  (see  Terzaghi  (1943))  requires  three 
nonparallel  subgroups  of  piles  in  the  foundation.  The  piles  within  each  subgroup 
are  assumed  to  be  parallel  and  are  assumed  to  have  the  same  head  load.  Each 
subgroup  is  replaced  by  a  single  pile  at  the  centroid  of  the  subgroup.  A  graphi¬ 
cal  procedure  is  used  to  resolve  the  superstructure  load  applied  to  the  rigid  cap  to 
each  subgroup. 


“Analytical”  Method 

Teng  (1962)  describes  a  simplified  procedure  for  including  the  effects  of 
horizontal  loads  as  well  as  battered  piles.  The  vertical  component  of  the  axial 
force  in  each  pile  due  to  the  resultant  vertical  load  and  moments  of  the  super¬ 
structure  on  the  rigid  cap  is  calculated  according  to  the  moment  of  inertia 
method.  The  total  axial  pile  load  and  its  horizontal  component  may  be  cal¬ 
culated  from  the  vertical  component.  Teng  suggests  that  an  adequate  design  has 
been  attained  if  the  applied  horizontal  foundation  load  does  not  exceed  the  sum 
of  the  horizontal  components  of  axial  pile  forces  by  more  than  1,000  lb/pile. 


Stiffness  Analysis  of  Pile  Foundations 

The  classical  methods  described  in  the  previous  paragraphs  essentially  neglect 
the  capability  of  the  piles  to  resist  lateral  loads  and  do  not  provide  a  means  of 
evaluating  the  stresses  induced  in  the  pile  by  bending  and  shear  at  the  pile  head. 
The  classical  methods  may  underestimate  the  strength  of  the  foundation  or  may 
lead  to  an  unconservative  design  depending  on  the  manner  in  which  the  pile 
head  is  attached  to  the  structure.  Hrennikoff  (1950)  and  Saul  (1968)  developed  a 
direct  stiffness  approach  to  the  analysis  of  two-  and  three-dimensional  pile 
groups  in  which  the  interaction  of  the  piles  with  the  surrounding  soil  as  well  as 
compatibility  of  pile  head  and  pile  cap  displacements  is  included.  In  this 
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procedure,  the  relationship  between  the  pile  head  forces  and  the  displacements 
of  the  point  of  attachment  to  the  rigid  pile  cap  is  assumed  to  be  linear. 


The  two  coordinate  systems  necessary  for  the  direct  stiffness  analysis  are  shown 
in  Figures  42  and  43  along  with  the  forces  and  displacements  on  the  pile  cap  and 
pile  head.  Relationships  between  the  global  and  local  axes  are  shown  in 
Figure  44. 


Figure  42.  Pile  cap  loads,  displacements,  and  coordinates 


The  pile  head  displacements  in  the  local  coordinate  system  for  a  pile  are 
expressed  in  terms  of  the  pile  cap  displacements  by  the  transformation 

{«}*  =  M,  [G\,  {U} 


where  {u}i  =  {«,  v;  w,  (pj  p;  0,}r  =  pile  head  displacements  in  the  local  coordinate 
system  for  the  z'th  pile;  [A\ ,  [G],  =  geometric  transformation  matrices  given  by 
(see  Figure  44  for  definitions  of  symbols) 
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Figure  43.  Head  forces,  displacements,  and  coordinates  for  iTH  pile 
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b.  View  A— A 


Figure  44.  Relationship  between  global  and  local  coordinates 
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and  {U}  =  {Uc  Vc  WG  (f)c  pG  dc}T  =  pile  cap  displacements  in  the  global  coordi¬ 
nate  directions. 

The  relationship  between  pile  cap  forces  and  pile  cap  displacements  is  given 
by 


{F} 


t  [G\-  [A]]  [b\  [A],  [G\t  J  { U } 


(146) 


where  {F}  =  {Fx  Fy  Fz  Mx  My  M,} T  =  superstructure  loads  on  pile  cap  in  global 
coordinate  directions;  n  =  number  of  piles  in  the  group;  and,  [6]j  =  pile  head 
stiffness  matrix  for  the  local  coordinate  system  described  in  Chapter  5. 
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After  Equation  146  has  been  solved  for  the  pile  cap  displacements,  the  head 
displacements  and  head  forces  in  the  local  coordinate  system  for  each  pile  are 
obtained  from  Equations  143  and  118,  respectively.  This  method  of  analysis  has 
been  incorporated  into  the  CASE  computer  program  CPGA  (Hartman  et  al. 

1984). 

The  pile-soil-pile  interaction  (PSPI)  method  is  identical  to  Saul's  method  in 
its  treatment  of  geometric  and  structural  analysis  aspects  of  the  pile  group 
problem.  However,  its  treatment  of  soil  reactions  against  piles  is  more  compre¬ 
hensive  and  rigorous.  The  method  was  first  formalized  by  O'Neill,  Ghazzaly  and 
Ha  (1977)  for  the  analysis  of  pile  clusters  in  offshore  drilling  platforms.  The 
PSPI  algorithm  PILGP2R  automates  the  process  of  computing  pile-head  stiff¬ 
ness  constants  bn  -  b66  and  bl5 ,  b24 ,  and  b42  by  first  developing  nonlinear 
relationships  between  axial  pile-head  load  and  axial  displacement,  pile-head 
torque  and  twist,  lateral  pile-head  load  and  displacement  with  zero  rotation,  pile- 
head  moment  and  rotation  with  zero  displacement,  and  the  cross-coupling 
relationship  between  rotation  and  lateral  load  and  between  displacement  and 
moment.  Sets  of  lateral  relationships  are  developed  for  each  of  two  orthogonal 
directions  to  accommodate  nonisotorepic  piles,  such  as  H-piles.  Nonlinear 
subroutines  similar  to  those  described  in  Chapters  2  and  3  are  used  to  develop 
the  axial  and  lateral  relationships,  respectively.  Hence,  the  user  inputs  become 
unit  load  transfer  curves.  Torsional  stiffness  of  individual  piles  is  generally  of 
minor  significance  in  large  groups  of  piles.  For  that  reason  special  computations 
of  torsional  stiffnesses  are  not  made,  and  the  torsional  stiffness  is  assumed  to  be 
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equal  to  a  constant  value  of  (GJ)pilJ2L  throughout  the  computations.  The  piles 
remain  linearly  elastic. 

Similarly, 

([/'),.  =  (cf  (A)  (148) 

where  (A)  is  the  vector  of  global  deflections  at  the  origin.  By  using  this  trans¬ 
formation,  each  pile  stiffness  is  moved  to  the  origin  of  the  pile  cap.  By  using 
Equations  93,  94,  and  95,  a  single-pile  force-deflection  relationship  can  be 


expressed  as 

(?),  =  (c),  (fl),  (a)f  ip)]  (A) 

(149) 

(0  =  E  i&i 

1  =  1 

(150) 

where 


(Q)  =  resultant  external  loads  applied  at  origin  of  pile  foundation 
n  =  total  number  of  piles  in  the  foundation 
This  leads  to  the  expression 


(0  =  h  w,  ifli  ipft  (a>  (i5i) 

/=1 

An  iterative  solution  is  used  to  determine  the  deformation  vector  for  the  rigid 
pile  cap.  Initial  b-factors  are  taken  as  the  initial  tangents  to  the  various  relation¬ 
ships  described  in  the  previous  paragraph.  The  three-dimensional  load  vector  is 
then  applied,  the  cap  deformation  vector  is  computed,  and  the  pile-head  defor¬ 
mation  vectors  are  computed.  From  these  deformations  pile-head  reactions  are 
computed  from  the  present  b-values.  A  check  is  then  made  of  the  compatibility 
of  the  reactions  and  deformations  at  each  pile  head  using  the  predetermined 
nonlinear  relationships.  If  compatibility  in  the  reaction  and  deformation  values 
is  noted,  new  b-values  are  computed  as  secants  to  the  axial  and  lateral  pile  head 
load-deformation  relationships  at  the  computed  values  of  deformation,  and  the 
process  is  repeated  as  necessary  until  tolerance  between  the  computed  deforma¬ 
tions  and  deformations  on  the  predetermined  relationships  at  the  computed  value 
of  pile-head  reaction  is  achieved  at  every  pile. 

Since  individual,  predetermined  nonlinear  relationships  are  used  to  represent 
the  relationships  between  the  various  lateral  deformation  modes  in  two  orthogo¬ 
nal  directions,  the  solution  that  is  thus  obtained  involves,  in  effect,  superposition 
of  nonlinear  functions,  which  is  not  theoretically  correct,  and  which  may 
produce  errors  of  practical  magnitude  as  the  pile-head  behavior  becomes  highly 
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nonlinear.  However,  where  nonlinearity  is  relatively  small,  as  would  typically 
occur  in  the  services  load  range,  the  errors  are  not  significant. 


Axial  and  lateral  behavior  algorithm  used 


Pile-soil-pile  interaction  is  then  included  as  indicated  in  Figure  45.  Once  the  . 
pile-head  deformation  is  known  from  the  preceding  solution,  the  axial  and  lateral 
behavior  algorithms  are  used  to  compute  the  soil  reaction  (axial  or  lateral)  at 
every  node  along  every  pile.  (In  practice,  only  selected  nodes  need  to  be  used  at 
this  stage  of  the  calculations.)  Lateral  reactions  are  computed  in  two  local 
orthogonal  directions.  The  lumped  reactions  against  the  soil  on  any  generic  Pile 
J  at  generic  Node  j  are  then  transformed  into  a  vector  of  orthogonal  reactions  in 
the  local  coordinate  system  (X,  Y,  Z)  shown  in  Figure  45.  The  additional 
displacements  in  the  local  X,  Y,  Z  directions  produced  by  these  reactions  on  any 
other  pile  1  in  the  system,  denoted  dace  ( k  =  X,  Y,  or  Z)  in  Figure  45,  are  then 
computed  by  using  Mindlin's  equations  for  displacements  in  the  interior  of  a 
semi-infinite  halfspace.  If  a  rigid  layer  (e.g.,  rock)  exists,  the  deflections 
computed  using  halfspace  theory  are  corrected  as  indicated  on  Figure  45. 
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Figure  45.  Geometric  definitions  for  computation 
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It  is  therefore  necessary  for  the  analyst  to  prescribe  values  of  Young's  modulus 
and  Poisson's  ratio  for  the  soil,  which  can  vary  linearly  from  the  soil  surface  to 
the  top  of  the  rigid  surface,  and  depth  to  the  rigid  surface.  As  with  the  elastic 
methods  for  axial  response  calculations,  these  moduli  should  generally  conform 
to  values  obtained  at  low  magnitudes  of  strain  (O'Neill  1983). 

The  additional  displacements  at  Node  i,  Pile  I,  due  to  reactions  at  all  nodes  on 
all  piles  other  than  Pile  I,  are  then  calculated  and  summed,  and  the  resulting 
added  deformations  at  Node  i,  Pile  I  due  to  PSPI  are  transformed  back  into  the 
normal  pile  coordinate  system.  These  displacements  are  then  compared  with  the 
displacements  computed  initially,  without  considering  PSPI.  Adjustments  are 
made  to  the  unit  load  transfer  curves  (f-w ,  Q-w  p-u  in  two  orthogonal  directions) 
to  force  the  unit  load  transfer  curves  to  pass  through  the  soil  displacement  value 
defined  by  the  sum  of  the  initial  displacement  <5  and  the  additional  displacement 
8'  in  the  appropriate  direction  at  the  value  of  soil  reaction  computed  initially. 
Specifically,  the  entire  curve  is  “stretched”  by  multiplying  all  displacement 
values  on  the  original  curve  input  by  the  analyst  by  (8 '  +  8)1 8.  This  process  is 
illustrated  in  Figure  46.  It  is  evident  that  in  geometrically  complex  groups  every 
unit  load  transfer  curve  on  every  pile  will  be  modified  in  a  different  manner.  It 
is  also  evident  that  this  procedure  does  not  reduce  the  maximum  unit  load 
transfer  that  can  occur  at  any  node  in  either  the  axial  or  lateral  direction;  thus,  it 
does  not  reduce  the  “capacity  efficiency”  of  the  group. 

The  solution  for  pile-cap  deformations  is  then  repeated  using  the  modified 
unit  load  transfer  curves,  which  produce,  in  general,  different  b  matrices  for  each 
pile  in  the  system  during  this  pass  through  the  solution.  The  regeneration  of  the 
load-deformation  relationships  for  this  second  pass  for  every  pile  in  the  system 
is  the  source  of  the  high  computational  effort  required  for  this  method  compared 
with  Saul's  method  and  which  makes  it  generally  impractical  for  the  design 
office.  Once  the  new,  compatible  pile  head  loads  and  deformations  have  been 
determined,  the  cap  deformation  vector  is  defined,  and  the  axial  and  lateral 
algorithms  are  entered  with  pile-head  deformations  as  boundary  conditions  to 
compute  moments,  shears  and  axial  thrusts  along  the  piles,  if  desired.  If  greater 
accuracy  is  desired,  the  correction  process  can  be  repeated;  however,  one  set  of 
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Unit  Load  Transfer 


Computed  reaction,  noninteractive  step 
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(input) 
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Figure  46.  Modification  of  unit  load  transfer  relationship  for  group  effects  at 
Node  i,  Pile  I 


corrections  of  the  unit-load  transfer  curves  appears  appropriate  for  most 
problems. 


Flexible  pile  cap  analysis 

The  pile  analysis  methods  described  above  assume  that  the  pile  cap,  or 
structure  base  slab,  is  rigid  in  comparison  to  the  stiffness  of  the  piles.  For  many 
structures,  such  as  U-frame  lock  monoliths,  this  is  not  a  valid  assumption,  and 
the  flexibility  of  the  base  slab  should  be  considered.  Currently  there  are  no 
special-purpose  programs  to  perform  this  type  analysis,  so  the  use  of  large 
general-purpose  programs  like  SAP  or  STRUDL  that  can  represent  a  flexible 
pile-cap  foundation  and  piles  must  be  used.  The  pile  element  used  in  the  rigid 
base  method  has  been  added  to  several  versions  of  the  SAP  and  STRUDL 
programs.  Flexible  base  analyses  have  already  been  performed  for  pile-founded 
structures  designed  by  the  Corps  of  Engineers.  A  more  detailed  report  on 
flexible  base  analysis  will  be  furnished  at  some  future  date. 


Nonlinear  analysis 

One  of  the  assumptions  made  in  the  rigid  base  analysis  method  is  that  a  pile 
can  be  represented  by  a  set  of  linear  stiffnesses.  The  actual  behavior  of  the  pile- 
soil  system  may  be  highly  nonlinear.  Some  existing  programs  are  capable  of 
nonlinear  analysis  of  a  structure  that  is  supported  by  only  a  few  piles.  One  such 
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program  which  may  be  able  to  economically  analyze  a  large  number  of  piles 
under  a  foundation  using  a  nonlinear  model  is  PILGP2R.  However,  for  large 
structures  supported  by  as  many  as  200  piles,  nonlinear  analysis  is  not  currently 
practical.  The  results  of  the  two  types  of  analyses  are  very  close.  In  O’Neil  and 
Tsao  (1984)  it  was  found  that,  while  nonlinearity  was  an  insignificant  effect,  a 
rational  procedure  must  be  applied  to  adjusting  the  subgrade  reaction  to  be  able 
to  account  for  PSPI  in  typical  Corps'  structures. 
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Appendix  A 

Linear  Approximations  for 
Load  Deformation  of  Axial 
Piles 


Linear  Elastic  Analyses 


For  a  linearly  elastic  representation  of  the  pile-soil  system,  the  governing 
differential  equation  is 

EA  —  -  Ez(z)w  =  0  (Al) 

d 2 

Z 

where  Ez(z)  is  the  stiffness  of  the  axial  load-transfer  mechanism.  A  typical  f-w 
curve  is  shown  in  Figure  Al.  At  any  displacement  w  the  nonlinear  f-w  curve 
may  be  replaced  by  a  linear  secant  kfz)  .  The  total  stiffness  of  the  load-transfer 
mechanism  appearing  in  Equation  Al  is  given  by 

Ez(z)  =2  Rkf  (z)  (A2) 


Soil  Stiffness  Constant  with  Depth 

Even  though  Equation  Al  is  representative  of  a  linearly  elastic  system,  a 
closed-form  solution  is  possible  only  if  the  side  friction  stiffness  is  constant  with 
depth.  While  this  condition  is  not  encountered  in  most  soil-pile  systems,  the 
solution  for  constant  stiffness  serves  to  indicate  characteristics  of  the  system 
behavior. 

For  E,  constant  with  depth,  the  solution  of  Equation  Al  may  be  written  as 
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Al 


Figure  A1 .  Typical  f-w  curve 


w 


p  t  ( 

Cl  cos  h  —  -  sin  h  — 
EA  [  Tz  Tz 


(A3) 


where 

P0  =  axial  compressive  force  at  pile  head 
Tz  =  EA/E 

Cl  =  constant  to  be  determined  from  conditions  at  the  pile  tip  (z=L) 
At  the  pile  tip,  three  alternative  conditions  are  presented: 
a.  Tip  reaction  equals  zero:  At  z  =  L  ,  Q  =  0  ,  hence 


cos  h  L/Tz  cos  h  Z ^ 
sin  h  L!TZ  sin  h  Zmax 


(A4) 


b.  Tip  displacement  equals  zero:  At  z  =  L  ,  w  =  0  ,  hence 


Cl 


sin  h  Z„m 

max 

cos  h  Z 

max 


(A5) 
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c.  Tip  elastically  restrained:  At  z  =  L ,  Q  =  K,w, ,  hence 

C1  =  C0S  *  +  r  Sin  h 

Sin  *  Zma*  +  r  C0S  Zm« 

where 

r  =  K/EzTz 

Kt  =  stiffness  of  elastic  resistance  at  pile  tip 
Zmax  ~  L/Tz 


(A6) 


Pile  Head  Axial  Stiffness 

Of  particular  interest  is  the  relationship  between  the  pile  head  force  P0  and 
axial  displacement  w0 .  For  convenience,  this  is  defined  by 


EA  wc 


where 

a.  Pile  tip  free 

=  C0S  *  Zma 

a°  sin  h 

b.  Pile  tip  fixed 

=  Sin  h_  Zmax 

a°  cos  h 

c.  Pile  tip  elastically  restrained 

=  C0S  h  Zmax 

0  sin  h  Zmax  +  r  C0S  k  Zma 


(A7) 


(A8) 


(A9) 


(A10) 


a0  is  plotted  for  various  values  of  Zmax  in  Figure  A2.  It  is  observed  that  for 
values  of  Zmax  greater  than  2,  conditions  at  the  tip  have  a  negligible  effect  on  the 
pile  head  force-displacement  relationship.  As  will  be  discussed  later,  the  pile 
head  axial  force-displacement  relationship  forms  a  part  of  the  elastic  pile  head 
stiffness  matrix  used  in  the  analysis  of  pile  groups.  The  axial  pile  head  stiffness 
coefficient  is  given  by 
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A3 


Figure  A2.  Axial  stiffness  coefficient  for  constant  soil  stiffness 


=  EA  _1_ 

Tz  a0 


(All) 


Nondimensional  Analysis 
for  Variable  Soil  Stiffness 

As  discussed  earlier,  the  stiffness  of  the  side  friction  displacement 
relationship  increases  with  depth.  For  axial  head  loads  in  the  working  load  range 
(e.g.,  one-half  of  the  ultimate  load),  it  has  been  found  that  the  equivalent  elastic 
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side  friction  increases  approximately  linearly  with  depth  for  normally 
consolidated  clays  and  as  the  square  root  of  the  depth  in  homogeneous  sands.  In 
general,  these  variations  may  be  expressed  as 

E2(z)  =  K0  +  Kfzn  (A12; 


where 

K0  =  elastic  stiffness  of  the  side  friction  effect  at  the  ground  surface 

Kf=  elastic  stiffness  coefficient  of  the  side  friction  effect  in  units  of  force  per 
unit  length  of  pile  per  unit  deflection 

n-  1  for  a  linear  variation  with  depth 

n  —  1/2  for  variation  as  the  square  root  of  depth 

Because  the  stiffness  of  the  side  friction  effect  is  a  function  of  the  strength  of  the 
soil,  K0  will  be  zero  for  sands.  Some  adhesion  of  clay  soil  may  occur  at  the 
ground  surface  and  K  for  clays  may  not  be  zero.  However,  it  is  likely  that 
installation  effects  will  minimize  adhesion  near  the  ground  surface  and  a 
conservative  estimate  is  obtained  for  K0  =  0  .  For  the  general  variation  with  K0  — 
0,  the  governing  differential  equation  is 

EA  -  Kf  z”  w  =  o  (A13) 

dz2 

Closed  form  solutions  of  Equation  A13  do  not  exist.  However,  nondimensional 
solutions  may  be  obtained  with  relatively  simple  numerical  techniques  as 
described  below. 

Following  the  procedures  described  by  Matlock  and  Reese  (1962),  the 
following  nondimensional  parameters  are  defined. 


z  =  — 

(A  14) 

Zmwc  =  ~ 

(A15) 

max  rj, 
z 

HZ)  -  M  wCz) 

(A  16) 

PoT z 

Kt 

K£TZTX 

(A  17) 
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A5 


(A  18) 


,—(71+20)  _  EA 
Z  "  Kf 

Substitution  of  the  nondimensional  parameters  into  the  governing  differential 
equation  yields 

-  Z"  a(Z)  =  0  (A  19) 

dZ 2 

Equation  A19  may  be  solved  for  a(Z)  by  any  numerical  technique  (e.g.,  finite 
differences).  From  the  solution  of  Equation  A 19,  the  relationship  between  pile 
head  force  and  head  displacement  is  obtained  as  a{Z= 0) ,  whence 


EA  wo 

a„  =  - 

Po  Tz 


(A20) 


and  the  pile  head  axial  stiffness  is 


Po 

Wo 


EA  _1_ 
T.  ao 


(A21) 


Values  of  the  stiffness  parameter  are  plotted  in  Figure  A3  for  a  linear 
variation  of  side  friction  stiffness  with  depth  and  in  Figure  25  for  side  friction 
stiffness  varying  as  the  square  root  of  depth.  These  figures  indicate  that  the  tip 
reaction  stiffness  has  negligible  effect  for  values  of  Zmax  greater  than  2. 


Evaluation  of  Side  Friction 
Stiffness  for  Piles  in  Sand 


Method  ESSF1 

Mosher  (1984)  suggested  a  secant  stiffness  computed  using  an  axial 
displacement  of  0. 1  in.  from  his  expression  for  the  f-w  curve.  The  resulting 
estimate  of  elastic  stiffness  for  side  friction  is 


Ez  = 


144  +  0.1^ 


(A22) 
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Figure  A3.  Axial  stiffness  coefficient  for  soil  stiffness  varying  linearly  with  depth 
where 

kf=  initial  slope  of  the  f-w  curve  given  in  Table  1  in  psf 
fmax  =  ultimate  side  friction  given  in  Figure  4  converted  to  units  of  psf 
C=  2 

R  =  circumference  of  the  pile 
The  resulting  units  of  E.  are  in  psi. 
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Figure  A4.  Axial  stiffness  coefficient  for  soil  stiffness  varying  as  square  root  of 
depth 

The  expression  for  E,  in  Equation  A22  includes  an  implicit  variation  with 
depth  embodied  in  the  curves  for  fmax  as  a  function  of  relative  depth  (see  Figure  4 
in  the  main  text).  An  examination  of  these  curves  suggests  that  the  variation  of 
E,  approximates  a  curve  proportional  to  the  square  root  of  depth.  To  utilize  the 
side  friction  stiffness  for  this  procedure,  the  actual  variation  of  E,  could  be  fitted 
to  a  curve  which  varies  as  the  square  root  of  depth  (e.g.  by  least  squares  fit). 
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Method  ESSF2 


The  procedure  of  Kraft,  Ray,  and  Cagawa  (1981)  for  f-w  curves  for  sand  may 
be  used  to  approximate  elastic  conditions  through  the  use  of  the  initial  slope  of 
the  f-w  curve.  From  Equation  6  in  the  main  text,  the  initial  slope  of  the  f-w  curve 
yields: 

a.  For  a  homogeneous  medium,  a  constant  side  friction  stiffness  equal  to 


Efz)  = 


GtC 


R  In 


2 L  (1  -  v) 


(A23) 


b.  For  a  medium  in  which  shear  modulus  varies  from  zero  at  the  ground 
surface  to  G,  at  the  pile  tip,  a  linearly  varying  side  friction  stiffness  equal 
to 


ESz) 


GtC 


2Z,(1  -  v) 


R  In 


z 

L 


R 


(A24) 


c.  For  a  soil  affected  by  pile  installation,  the  shear  modulus  G  or  Gt 

appearing  in  Equations  A23  and  A24  should  be  replaced  by  an  effective 
modulus  as  given  in  Equation  10  or  1 1  in  the  main  text. 


Evaluation  of  Side  Friction 
Stiffness  for  Piles  in  Clay 

Method  ECSF1 

Kraft,  Ray,  and  Kagawa  (1981)  suggest  that  the  same  procedures  used  for 
sand  may  be  used  for  piles  in  clay  with  the  shear  moduli  appearing  in 
Equations  A23  and  A24  being  modified  to  account  for  magnitude  and  time  of 
loading,  as  shown  in  Table  Al. 


Method  ECSF2 

Heydinger  (1984)  suggests  using  a  secant  to  the  f-w  curve  evaluated  at  a 
relative  displacement  w/2R  =  0.005.  From  Equation  22  in  the  main  text, 
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Table  A1 

Adjustment  in  Gfor  Various  Loading  Conditions  (Adjustment 
factor  =  G  (operational/G  (in  situ)) 

Condition 

Adjustment  Factor 

Immediately  After  Driving 

Small  load 

0.85 

Large  load  (>  working  load) 

0.54 

After  Setup 

Small  load 

1.04 

Large  load  (f  =  0Asu) 

0.96 

Large  load  (f  =  $,) 

0.73 

Ez- 


Ef  C/2R 


1  + 


0.005- 


Ef 


fr 


max  ) 


(A25) 


where  island  m  are  given  by  Equations  23  and  24  in  the  main  text,  respectively. 
The  variation  of  E2  along  the  pile  depends  on  the  distribution  of  soil  modulus  of 
elasticity  used  in  the  evaluation  of  Eyfrorn  Equation  23. 


Method  ECSF2 


The  elasto-plastic  representation  of  side  friction  due  to  Aschenbrener  and 
Olson  (1984)  yields 


Ez  =  — - - 

z  0.1  (in.) 

where  for  consistent  units  su  should  be  expressed  in  psi  and  C  in  inches 
producing  E,  in  psi. 


(A26) 


Evaluation  of  Tip  Reaction  Stiffness 

General 

As  shown  previously,  the  tip  reaction  only  has  a  significant  effect  on  the  pile 
head  stiffness  coefficient  for  piles  having  Zmax  less  than  2.  If  the  value  of  Zmm 
resulting  from  any  of  the  assessments  of  side  friction  described  above  is  less 
than  2,  the  tip  reaction  stiffness  may  be  omitted. 
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In  the  discussions  for  evaluating  the  tip  reaction  stiffness  that  follow, 
stiffness  is  proportional  to  the  effective  area  at  the  tip  At  bearing  on  the  soil.  For 
closed  end  or  solid  piles  the  effective  tip  area  may  reasonably  be  taken  as  the 
cross  section  area  of  the  pile.  For  H-piles  or  open-ended  pipe  piles  the  tip  area 
may  be  as  little  as  the  area  of  material  in  the  cross  section  to  an  area  equal  to  that 
bounded  by  the  exterior  of  the  section  (see  Figure  5  in  the  main  text).  When  the 
radius  of  the  tip  reaction  area  is  required  to  evaluate  tip  stiffness,  an  effective 
radius  is  obtained  from 


(A27) 


The  tip  reaction  stiffness  may  be  obtained  from  any  of  the  procedures 
described  previously  for  developing  q-w  curves  by  evaluating  a  secant  stiffness 
for  a  tip  displacement  representative  of  working  load  conditions.  Typically  in 
the  Corps  of  Engineers,  failure  at  the  tip  is  considered  to  occur  at  a  tip 
displacement  of  0.25  in.  Unless  stated  otherwise,  working  load  conditions  are 
assumed  to  occur  at  one-tenth  of  the  displacement  corresponding  to  failure  (i.e., 
0.025  in.). 


Evaluation  of  Tip  Reaction 
Stiffness  for  Piles  in  Sand 


Method  EST1 

The  theory  of  elasticity  solution  for  a  rigid  punch  has  been  used  by  Kraft, 
Ray,  and  Kagawa  (1981)  (see  also  Randolph  and  Wroth  1978)  to  estimate  the  tip 
reaction  stiffness  as 


Kt  = 


ARtG 
7,(1  -  v) 


(A28) 


where  the  shear  modulus  G  should  be  taken  as  an  average  in  situ  value  between 
6 Rt  above  the  pile  tip  to  6 Rt  below  the  tip.  The  factor  /,  in  Equation  A28  is  an 
influence  factor  ranging  from  0.5  to  0.78. 


Method  EST2 

Mosher  (1984)  and  Vijayvergiya  (1977)  express  the  tip  reaction  q-w  curve  as 
a  power  function  (see  pages  24-26).  Mosher  recommends  for  working  load 
approximations  a  secant  tip  reaction  stiffness  corresponding  to  a  tip 
displacement  of  0.025  in.  The  corresponding  tip  stiffnesses  are: 

a.  For  loose  sand:  K,  =  12.6  A ,  qmax  (A29) 

b.  For  medium  sand:  K,  =  18.6  At  qmax  (A30) 
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c.  For  dense  sand:  K,  -  22.5  A,  qmax  (A3 1) 

where  qmca  is  the  ultimate  unit  tip  reaction  from  Figure  17.  For  consistent  units 
in  Equations  A29  through  A3 1,  qmax  must  be  in  pounds  per  square  inch,  and  A, 
must  be  in  square  inches,  which  yields  Kt  in  pounds  per  inch. 


Method  EST3 

A  secant  stiffness  obtained  from  the  work  of  Briaud  and  Tucker  (1984), 
which  considers  the  effects  of  residual  stresses  due  to  installation  for  a  tip 
displacement  0.0.025  in.,  is 


g  _  KflmaAl 

'  +  0.025*, 

"max  q 

and 

kq  =  467.1  TV0  0065 


(A32) 


(A33) 


where  N  is  the  average  uncorrected  standard  penetration  count  in  blows  per  foot 
from  a  distance  of  87?,  above  the  pile  tip  to  8 7?,  below  the  tip.  The  units  of  kq  in 
Equation  A33  are  tsf/in.  The  required  units  of  other  terms  in  Equation  A32  are 
qmax  in  tons  per  square  foot  and  A,  in  square  feet,  which  yields  Kt  in  tons  per  inch. 


Evaluation  of  Tip  Reaction 
Stiffness  for  Piles  in  Clay 

Method  ECT1 

The  bilinear  tip  reaction  curve  used  by  Aschenbrener  and  Olson  (1984) 
produces 


=  9lA 

0.027?, 


(A34) 


where  su  is  the  average  undrained  shear  strength  of  the  clay  from  67?,  above  the 
pile  tip  to  67?,  below  the  tip. 


Method  ECT2 

The  tip  stiffness  developed  by  Kraft,  Ray,  and  Kagawa  (1981)  described  on 
page  All  may  be  used  for  piles  in  clay. 
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Method  ECT3 


Skempton  (1951)  observed  the  similarity  of  the  load-displacement  behavior 
of  a  plate  load  test  and  the  laboratory  stress-strain  curve  for  soft  clays.  It  was 
concluded  that  a  linear  approximation  of  the  load  displacement  relationship  up  to 
half  of  the  ultimate  load  could  be  related  to  the  strain  at  50  percent  of  the 
unconfmed  compression  strength  indicated  by  the  laboratory  stress-strain  curve. 
The  observation  has  been  used  to  obtain  an  estimate  of  the  pile  tip  reaction 
stiffness  as 


Kt  = 


^e50 


(A35) 


where 

qu  =  unconfmed  compression  strength  of  the  clay  at  the  pile  tip 
A,  =  effective  tip  area 

e50  =  strain  at  50  percent  of  ultimate  strength  from  a  laboratory  stress-strain 
curve 

R,  =  effective  radius  of  the  tip  area 

Typical  values  of  e50  are  0.02  for  a  very  soft  clay,  0.01  for  a  soft  clay,  and  0.005 
for  a  stiff  clay. 
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Appendix  B 

Nondimensional  Coefficients 
for  Laterally  Loaded  Piles 


Basic  Equations 

Ez  =  kzn 

jyn+4  _  EI 

k 

EI  dMz)_  +  =  Q 

dz 4 


L 

T 


u{Z) 


VT3 

-JTa“(Z) 


M0T2 

EI 


BU(Z) 


du 

dz 


VT2  MJ 

—  + 


EI 


—  B(Z) 
EI  5 


M(Z)  =  VJAJZ)  +  M0BJZ) 
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Table  B1 

Nondimensional  Coefficients  for  Laterally  Loaded  Pile  for  Soil  Modulus 
Depth  (Head  Shear  Va  =  1 ,  Head  Moment  Mn  =  0 


Constant  with 


1.597 

-1.116 

0.000 

1.474 

-1.032 

0.000 

1.486 

-1.111 

0.094 

1.371 

-1.027 

0.095 

1.375 

-1.097 

0.174 

1.248 

-1.009 

0.192 

1.266 

-1.076 

0.240 

1.168 

-0.992 

0.247 

1.160 

-1.050 

0.293 

1.051 

-0.958 

0.315 

1.057 

-1.018 

0.335 

0.976 

-0.931 

0.351 

0.956 

-0.983 

0.366 

0.867 

-0.886 

0.395 

0.860 

-0.945 

0.388 

0.797 

-0.854 

0.417 

0.767 

-0.906 

0.401 

0.698 

-0.802 

0.441 

0.679 

-0.865 

0.406 

0.635 

-0.767 

0.451 

0.594 

-0.825 

0.405 

0.546 

-0.712 

0.458 

0.514 

-0.784 

0.397 

0.491 

-0.675 

0.459 

0.437 

-0.745 

0.385 

0.413 

-0.620 

0.454 

0.365 

-0.708 

0.368 

0.365 

-0.584 

0.447 

0.296 

-0.672 

0.347 

0.298 

-0.532 

0.433 

0.230 

-0.638 

0.324 

0.257 

-0.497 

0.421 

0.168 

-0.607 

0.298 

0.200 

-0.448 

0.400 

0.109 

-0.579 

0.270 

0.165 

-0.417 

0.385 

0.052 

-0.553 

0.241 

0.118 

-0.372 

0.359 

-0.002 

-0.531 

0.212 

0.089 

-0.344 

0.341 

-0.054 

-0.511 

0.183 

0.051 

-0.305 

0.314 

-0.104 

-0.494 

0.154 

0.027 

-0.280 

0.295 

-0.153 

-0.480 

0.127 

-0.004 

-0.247 

0.266 

-0.200 

-0.469 

0.101 

-0.023 

-0.226 

0.247 

-0.247 

-0.460 

0.077 

-0.049 

-0.198 

0.218 

-0.292 

-0.453 

0.055 

-0.064 

-0.182 

0.200 

-0.337 

-0.449 

0.037 

-0.084 

-0.159 

0.173 

-0.382 

-0.446 

0.021 

-0.097 

-0.146 

0.155 
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Table  B1  (Concluded) 

z 

=  2 

2  max  ~  ® 

2^=4  | 

A  A 

Mu  ....  "s _ 

Am  . 

n  ■ 

As 

4* 

Au 

A  A 

2.80 

-0.426 

-0.444 

0.010 

-0.113 

-0.129 

0.130 

2.90 

-0.471 

-0.444 

0.003 

-0.123 

-0.119 

0.115 

3.00 

-0.515 

-0.443 

0.000 

-0.136 

-0.107 

0.093 

3.10 

-0.145 

-0.100 

0.079 

3.20 

-0.156 

-0.091 

0.060 

3.30 

-0.163 

-0.087 

0.049 

3.40 

-0.173 

-0.082 

0.034 

3.50 

-0.180 

-0.080 

0.026 

3.60 

-0.189 

-0.077 

0.015 

3.70 

-0.195 

-0.076 

0.010 

3.80 

-0.205 

-0.076 

0.004 

3.90 

-0.21 1 

-0.075 

0.001 

4.00 

-0.218 

-0.075 

0.000 
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Table  B2 

Nondimensional  Coefficients  for  Laterally  Loaded  Pile  for  Soil 
Depth  (Head  Shear  Va  =  0,  Head  Moment  Mn=1) _ 


Modulus  Constant  with 


-1.741 


-1.519 


-1.308 


-1.107 


-0.915 


-0.731 


-0.557 


-0.389 


-0.229 


-0.076 


0.072 


0.216 


0.355 


0.491 


0.624 


0.755 


0.885 


1.013 


1.141 


1.269 


1.397 


2.261 


2.161 


2.063 


1.967 


1.876 


1.789 


1.709 


1.635 


1.568 


1.508 


1.456 


1.412 


1.374 


1.344 


1.320 


1.303 


1.290 


1.283 


1.279 


1.277 


1.277 


-1.000 


-0.993 


-0.972 


-0.937 


-0.891 


-0.836 


-0.773 


-0.705 


-0.633 


-0.559 


-0.484 


-0.410 


-0.338 


-0.269 


-0.205 


-0.148 


-0.098 


-0.057 


-0.026 


-0.007 


0.000 


-1.116 

1.500 

-1.000 

-1.032 

1.441 

-0.971 

1.400 

-0.996 

-0.893 

1.341 

-0.997 

-0.836 

1.301 

-0.983 

-0.739 

1.222 

-0.982 

-0.710 

1.204 

-0.961 

-0.645 

1.144 

-0.966 

-0.595 

1.109 

-0.932 

-0.514 

1.030 

-0.934 

-0.489 

1.017 

-0.897 

-0.435 

0.956 

-0.909 

-0.391 

0.930 

-0.858 

-0.327 

0.850 

-0.865 

-0.303 

0.846 

-0.814 

-0.261 

0.782 

-0.834 

-0.222 

0.767 

-0.767 

-0.173 

0.685 

-0.784 

-0.149 

0.693 

-0.719 

-0.121 

0.623 

-0.749 

-0.083 

0.623 

-0.668 

-0.052 

0.537 

-0.695 

-0.024 

0.559 

-0.617 

-0.01 1 

0.483 

-0.658 

0.029 

0.500 

-0.565 

0.042 

0.407 

-0.604 

0.076 

0.446 

-0.514 

0.073 

0.360 

-0.568 

0.118 

0.397 

-0.464 

0.112 

0.295 

-0.515 

0.156 

0.353 

-0.415 

0.134 

0.255 

-0.480 

0.189 

0.314 

-0.367 

0.162 

0.201 

-0.430 

0.218 

0.280 

-0.321 

0.176 

0.168 

-0.397 

0.245 

0.250 

-0.277 

0.194 

0.123 

-0.351 

0.269 

0.224 

-0.236 

0.202 

0.096 

-0.322 

0.290 

0.202 

-0.198 

0.212 

0.060 

-0.280 

0.309 

0.184 

-0.162 

0.216 

0.038 

-0.254 

0.327 

0.170 

-0.130 

0.218 

0.010 

-0.218 

0.343 

0.158 

-0.100 

0.219 

-0.006 

-0.196 

0.359 

0.150 

-0.075 

0.216 

-0.028 

-0.164 

0.373 

0.143 

-0.052 

0.214 

-0.040 

-0.145 

0.387 

0.139 

-0.034 

0.208 

-0.056 

-0.119 

0.401 

0.136 

-0.019 

0.203 

-0.065 

-0.104 

( Continued 
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1  Table  B2  (Concluded) 

2^=2 

=  3 

Zm.r=4  I 

z 

Bs 

Bu  

** 

Bm  . 

Bu 

B* 

*0, 

2.80 

0.415 

0.135 

-0.009 

0.195 

-0.076 

-0.083 

2.90 

0.428 

0.134 

-0.002 

0.188 

'  -0.082 

-0.070 

3.00 

0.442 

0.134 

0.000 

0.178 

-0.090 

-0.054 

3.10 

0.170 

-0.094 

-0.044 

3.20 

0.159 

-0.098 

-0.032 

3.30 

0.151 

-0.101 

-0.025 

3.40 

0.139 

-0.103 

-0.017 

3.50 

0.130 

-0.104 

-0.012 

3.60 

0.118 

-0.105 

-0.007 

3.70 

0.109 

-0.106 

-0.004 

3.80 

0.097 

-0.106 

-0.001 

3.90 

0.088 

-0.106 

-0.000 

4.00 

0.078 

-0.106 

0.000 
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Figure  B4.  Shear  coefficient  for  unit  head  shear  for  soil  stiffness  constant  with  depth 
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Figure  B6.  Slope  coefficient  for  unit  head  shear  for  soil  stiffness  varying  linearly  with  depth 
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Table  B3 

Nondimensional  Coefficients  for 
Linearly  with  Depth  (Head  Shear 


Laterally  Loaded  Pile  for  Soil  Modulus  Varying 
Va  =  1 ,  Head  Moment  Ma  =  0) 


2.442 

-1.622 

0.000 

2.280 

-1.616 

0.100 

2.087 

-1.597 

0.216 

1.960 

-1.577 

0.290 

1.773 

-1.536 

0.394 

1.651 

-1.502 

0.458 

1.474 

-1.442 

0.543 

1.361 

-1.396 

0.592 

1.198 

-1.321 

0.655 

1.094 

-1.268 

0.689 

0.947 

-1.182 

0.729 

0.855 

-1.123 

0.747 

0.725 

-1.032 

0.764 

0.645 

-0.971 

0.768 

0.534 

-0.879 

0.764 

0.466 

-0.818 

0.755 

0.374 

-0.729 

0.733 

0.318 

-0.671 

0.714 

0.242 

-0.587 

0.679 

0.197 

-0.534 

0.652 

0.138 

-0.459 

0.607 

0.103 

-0.411 

0.574 

0.058 

-0.346 

0.523 

0.032 

-0.305 

0.488 

-0.002 

-0.250 

0.435 

-0.020 

-0.216 

0.399 

-0.043 

-0.172 

0.346 

-0.056 

-0.145 

0.312 
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Table  B3  (Concluded) 

z 

S-r«  3 

=  4  II 

A„ 

A . 

A 

A. 

4. 

* 

Am   

2.80 

-0.383 

-0.559 

0.027 

-0.071 

-0.111 

0.262 

2.90 

-0.439 

-0.558 

0.007 

-0.079 

-0.091  * 

0.231 

3.00 

-0.494 

-0.557 

0.000 

-0.089 

-0.066 

0.186 

3.10 

-0.093 

-0.052 

0.159 

3.20 

-0.099 

-0.036 

0.121 

3.30 

-0.101 

-0.027 

0.098 

3.40 

-0.104 

-0.017 

0.069 

3.50 

-0.105 

-0.012 

0.051 

3.60 

-0.106 

-0.007 

0.030 

3.70 

-0.106 

-0.005 

0.019 

3.80 

-0.107 

-0.004 

0.007 

3.90 

-0.107 

-0.003 

0.002 

4.00 

-0.108 

-0.003 

0.000 
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Figure  B8.  Shear  coefficient  for  unit  head  shear  for  soil  stiffness  varying  linearly  with  depth 
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Figure  B9.  Deflection  coefficient  for  unit  head  shear  for  soil  stiffness  varying  linearly  with  depth 
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Table  B4 

Nondimensional  Coefficients  for 
Linearly  with  Depth  (Head  Shear 


Laterally  Loaded  Pile  for  Soil  Modulus  Varying 
VD  =  0,  Head  Moment  =  1) _  _ 
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Table  B4  (Concluded) 

=  2 

Zm.,=  3 

_ _ II 

z 

B. . 

Bs 

Bm  .... 

Bu 

Bs 

B.  ...  . . - 

Bu 

Bs 

Bm 

2.80 

0.259 

0.163 

-0.017 

0.105 

-0.046 

-0.115 

2.90 

0.276 

0.162 

-0.004 

0.101 

-0.054 

-0.097 

3.00 

0.292 

0.162 

0.000 

0.093 

-0.064 

-0.072 

3.10 

0.088 

-0.069 

-0.059 

3.20 

0.079 

-0.075 

-0.041 

3.30 

0.073 

-0.078 

-0.031 

3.40 

0.064 

-0.081 

-0.020 

3.50 

0.057 

-0.083 

-0.014 

3.60 

0.047 

-0.084 

-0.007 

3.70 

0.040 

-0.084 

-0.004 

3.80 

0.030 

-0.085 

-0.001 

3.90 

0.023 

-0.085 

-0.000 

4.00 

0.015 

-0.085 

-0.000 
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Relative  Depth 


Figure  B1 1 .  Bending  moment  coefficient  for  unit  head  shear  for  soil  stiffness  varying 
parabolically  with  depth 
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Figure  B1 2.  Shear  coefficient  for  unit  head  shear  for  soil  stiffness  varying  parabolically  with 
deDth 
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Table  B5 

Nondimensional  Coefficients  for  Laterally 
Parabolically  with  Depth  (Head  Shear  V.  = 


Loaded  Pile  for  Soil  Modulus  Varying 
1,  Head  Moment  Ma  =  0) 
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Table  B5  (Concluded) 

Z  =2 

*-max  *- 

Z  =3 

•-max  w 

2™,  =  4 

z 

A„ 

As 

Am 

H 

_Am - 

A. 

A. 

A. 

2.80 

-0.221 

-0.394 

0.047 

-0.052 

-0.050 

0.298 

2.90 

-0.260 

-0.391 

0.013 

-0.055 

-0.028 

0.256 

3.00 

-0.299 

-0.390 

0.000 

-0.057 

-0.001 

0.198 

3.10 

-0.057 

0.013 

0.163 

3.20 

-0.054 

0.030 

0.117 

3.30 

-0.051 

0.038 

0.091 

3.40 

-0.046 

0.047 

0.059 

3.50 

-0.042 

0.051 

0.042 

3.60 

-0.036 

0.055 

0.023 

3.70 

-0.031 

0.057 

0.013 

3.80 

-0.024 

0.058 

0.004 

3.90 

-0.020 

0.058 

0.001 

4.00 

-0.014 

0.058 

0.000 
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Figure  B13.  Deflection  coefficient  for  unit  head  moment  for  soil  stiffness  constant  with  depth 
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Figure  B14.  Slope  coefficient  for  unit  head  moment  for  soil  stiffness  constant  with  depth 
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Figure  B1 5.  Bending  moment  coefficient  for  unit  head  moment  for  soil  stiffness  constant  with 
depth 
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Table  B6 

Nondimensional  Coefficients  for  Laterally 
Parabolical  I  v  with  Depth  (Head  Shear  Vn  = 


Loaded  Pile  for  Soil 
0,  Head  Moment 


Modulus  Varying 
=  1) 
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Table  B6  (Concluded) 

z^  =  2 

3 

_ _ II 

z 

Bu 

Bs 

Bu 

e. 

B 

B 

B. . 

.B»~  - 

2.80 

0.126 

0.093 

-0.024 

0.056 

-0.035 

-0.103 

2.90 

0.136  * 

0.092 

-0.006 

0.053 

-0.042 

-0.081 

3.00 

0.145 

0.092 

0.000 

0.048 

-0.050 

-0.053 

3.10 

0.043 

-0.054 

-0.038 

3.20 

0.037 

-0.057 

-0.021 

3.30 

0.032 

-0.059 

-0.012 

3.40 

0.025 

-0.060 

-0.003 

3.50 

0.020 

-0.060 

0.000 

3.60 

0.013 

-0.060 

0.002  ; 

3.70 

0.008 

-0.059 

0.003 

3.80 

0.001 

-0.059 

0.002 

3.90 

-0.004 

-0.059 

0.001 

4.00 

-0.009 

-0.059 

0.000 
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Figure  B1 6.  Shear  coefficient  for  unit  head  moment  for  soil  stiffness  constant  with  depth 
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Figure  B1 9.  Bending  moment  coefficient  for  unit  head  moment  for  soil  stiffness  varying  linearly 
with  depth 
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Figure  B20.  Shear  coefficient  for  unit  head  moment  for  soil  stiffness  varying  linearly  with  depth 
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Figure  B23.  Bending  moment  coefficient  for  unit  head  moment  for  soil  stiffness  varying 
parabolically  with  depth 
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Figure  B25.  Pile  head  deflection  coefficients  for  unit  head  shear 
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Figure  B26.  Pile  head  slope  coefficients  for  unit  head  shear 
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Figure  B28.  Pile  head  slope  coefficients  for  unit  head  moment 
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